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1.  INTRODUCTION 

1.1.  General 

This  Final  Report  describes  the  complete  research  work  that  has  been  carried  out  starting 
1  April  1995,  under  Office  of  Naval  Research  SBIR  Phase  I  contract  no.  N00014-95-C-0154,  in 
the  area  of  energetics  of  late  chemical  reactions  of  metals  in  nonideal  underwater  detonations. 

1.2.  Bacl^round 

The  reaction  of  metallic  fuels  with  the  detonation  products  in  all  explosives  is  poorly 
understood  to  this  day,  but  there  is  much  evidence,  albeit  indirect,  that  the  exothermic  reaction 
that  consmnes  aluminum  takes  place  over  a  much  longer  time  scale  than  that  which  consumes, 
say,  TNT,  RDX,  or  HMX.  The  probability  of  this  late  time  reaction  is  greater  in  underwater 
explosions  because  the  inertial  confinement  of  the  explosive  products  by  the  vast  mass  of 
surrounding  water  is  much  greater  than  any  confinement  by  metal  cases  of  explosions  in  air. 
Metallic  fuels,  principally  aluminum,  have  been  added  to  underwater  explosives,  and  the 
aluminum  oxidation  is  greatly  affected  by  the  confinement.  It  is  well  known  that  the  "bubble 
energy"  -  a  uniquely  underwater  effect  ~  increases,  relative  to  the  nonaluminized  material,  as  the 
aluminum/oxygen  ratio  increases  from  0  to  1 .  The  "shock  wave  energy"  reaches  a  maximum  for 
an  Al/0  ratio  of  about  0.4  in  some  TNT/RDX  mixtures.  The  aluminum  content  is  used  to  fine 
tune  the  bubble  and  shock  energies  -  which  are  related  to  underwater  target  damage  -  depending 
on  the  target. 

Metallized  composite  explosives  provide  the  means  of  tailoring  the  chemical  energy 
release  rate  in  the  detonation  process  because  the  metal  additive  is  capable  of  reacting  with  the 
explosive's  products  to  support  both  the  propagation  of  the  detonation  wave  and  the  unsteady 
flow  in  the  release  wave  behind  the  detonation  front.  Because  of  this  role  of  the  metal  in  the 
detonation  process,  metallized  explosives  are  said  to  be  nonideal  explosives;  a  consequence  is 
that  there  are  not  predictive  models  for  quantifying  their  behavior.  In  addition,  because  the 
amount  of  metal  reacting  to  support  the  detonation  is  not  known,  and  thermo-hydrodynamic 
codes  for  calculating  Chapman-Jouguet  (CJ)  parameters  for  CHNO  explosives  cannot  be  used  to 
predict  detonation  parameters  in  metallized  explosives,  the  detonation  process  is  also  said  to  be 
nonideal.  Consequently,  the  capability  to  tailor  the  energy  release  in  metallized  explosives  for 
underwater  applications  depends  on  our  ability  to  control  the  reactions  of  the  metal  in  the 
nonideal  detonation  process. 

Metallized  explosives  are  examples  of  composite  explosives  that  are  nonideal  explosives, 
i.e,  where  mass  transport  between  oxidizer  and  fuel  is  an  important  rate-determining  process, 
unlike  ideal  explosives,  which  are  molecular  mixtures  of  fuel  and  oxidizer  that  are  separated  by 
extremely  short  distances,  i.e.,  bond  lengths.  Because  reaction  rates  in  ideal  explosives  are  so 
fast,  explosive  performance  can  be  modeled  as  depending  only  on  equilibrium  states  along  the 
Chapman-Jouguet  (CJ)  isentrope  behind  the  CJ  plane,  upstream  of  which  all  explosive  is  reacted. 
The  modeling  of  nonideal  detonations,  i.e.,  detonations  of  nonideal  explosives,  is  much  more 
complicated,  depending  on  the  heterogeneity  of  the  fuel/oxidizer  mixture,  diffusion  rates,  mixing 
distances,  reaction  times,  and  energy  release  rates,  none  of  which  are  explicitly  needed  to 
describe  ideal  detonations. 
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Ideal  and  nonideal  explosives  differ  greatly  in  their  failure  diameter:  failure  diameters  in 
the  former  are  generally  of  the  order  of  millimeters;  those  of  the  latter  are  in  the  centimeter  range 
because  of  the  heterogeneity  and  slow  reaction  rates  of  the  composite  explosives. 

The  energy  release  processes  in  composites,  while  relatively  slow,  are  far  more  complex 
than  in  ideal  explosives.  As  a  consequence,  it  is  a  difficult  problem  to  couple  these 
"microscopic"  rate-determining  processes  to  the  "macroscopic"  compressible  fluid  dynamics  of 
the  complex  detonation  products  and  the  "early"  compressible  and  "later"  incompressible  flows 
in  the  water. 

For  ideal  and  nonideal  detonation  processes,  the  flows  in  the  steady-state  reaction  zones 
satisfy  the  Rankine-Hugoniot  (RH)  jump  conditions  and  become  sonic  at  the  CJ  point  where  the 
energy  release  rate  becomes  zero.'-^  But  whereas  all  the  available  chemical  energy  is  used  to 
support  the  steady-state  flow  in  the  ideal  detonation,  it  is  not  used  in  the  nonideal  detonation 
because  the  time  scale  for  the  decomposition  of  the  explosive  is  significantly  different  firom  the 
time  scales  of  the  reactions  of  the  metal  with  the  detonation  products.  The  detonation  parameters 
in  the  nonideal  detonation  thus  can  be  tailored  by  controlling  the  amount  of  metal  that  reacts  in 
the  steady-state  reaction  zone. 

For  both  detonation  processes,  the  flows  in  the  unsteady  release  waves  behind  the  CJ 
planes  are  governed  by  the  differential  equations  expressing  the  conservation  of  mass, 
momentum,  and  energy  for  inviscid,  adiabatic  flow.  The  release  wave  in  the  ideal  detonation 
process  is  isentropic  and  is,  thus,  represented  in  the  specific  volume-pressure  v-p  plane  by  the 
isentrope  of  the  detonation  products  passing  through  the  CJ  point.  Because  of  metal  reactions, 
however,  the  release  wave  in  the  nonideal  detonation  process  is  not  isentropic  and  is,  thus, 
represented  in  the  v-p  plane  by  a  reactive  release  adiabat  passing  through  the  CJ  point. 

Moreover,  the  entropy  increases  with  volume  along  such  a  reactive  adiabat  because  entropy  is 
produced  by  the  reactions  of  the  metal  with  the  products  of  the  explosive.  Thus,  the  energy 
release  rate  for  underwater  applications  can  be  tailored  by  controlling  the  reactions  of  the  metal 
with  the  explosive's  products  in  the  release  wave. 

Because  the  areas  imder  the  CJ  isentrope  and  CJ  release  adiabat  in  the  v-p  plane  are 
directly  related  to  the  amounts  of  available  chemical  energy  in  the  explosive  compositions,  the 
energy  release  rate  in  the  release  wave  can  be  tailored  by  controlling  the  pressure  when  the  metal 
begins  to  react  as  well  as  its  subsequent  rate  of  reaction.  In  other  words,  the  shapes  of  the 
reactive  release  adiabats  can  only  be  changed  by  varying  both  the  onset  and  subsequent  rates  of 
the  metal  reactions  in  the  release  wave. 

It  follows  from  the  discussion  of  the  detonation  process  that  ideal  detonations  can  be 
adequately  described  by  thermodynamic  processes,  but  nonideal  detonation  caimot  because  the 
energy  release  rate  is  governed  by  kinetic  processes  with  different  time  scales.  Semi-empirical 
thermo-hydrodynamic  studies  of  nonideal  explosives,  of  the  type  performed  by  Mader^  and 
Cowperthwaite,'*  provide  phenomenological  information  about  the  energy  release  rate,  but  do  not 
provide  the  information  and  knowledge  required  to  control  it.  Thus,  a  more  basic  approach  must 
be  taken  to  provide  a  means  of  controlling  the  energy  release  rate  in  metallized  explosives.  In 
such  an  approach,  we  must  first  identify  the  most  important  kinetic  processes  governing  the 
energy  release  rate  and  then  develop  realistic  models  for  them,  containing  parameters  related  to 
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the  physical  and  chemical  properties  of  the  explosive  compositions. 

1.3.  Dividing  Up  the  Kinds  of  Explosives 

Explosives  are  generally  divided  into  two  classes,  ideal  and  nonideal.^  We  have  noted 
that  ideal  explosives  may  be  defined  as  those  in  which  fuel  and  oxidizer  are  homogeneously 
mixed  at  the  molecular  level.  That  is  to  say,  the  irreducible  fuel  element  and  the  irreducible 
oxidizer  element  are  either,  within  the  same  molecule  in  a  grain  or  crystal  of  identical  molecules, 
or,  are  at  least  nearest  neighbors  in  separate  molecules  in  a  bimolecular  crystal  or  grain.  Ideal 
explosives  are  characterized  by  short  reaction  zones  and  exhibit  little  or  no  delayed  reaction. 
Nonideal  explosives  are  a  broader  class  usually  composed  of  grains  of  ideal  explosive 
homogeneously  mixed  with  other  grains  of  fuel  and/or  oxidizer.  They  may  behave  more  or  less 
like  ideal  explosives  depending  on  whether  they  contain  more  or  less  ideal  explosive 
constituents.  They  are  composite  explosives,  exhibit  longer  reaction  zones,  and  in  general 
considerably  more  "delayed"  reaction. 

For  the  purpose  of  modeling  nonideal  explosives,  we  divide  these  into  two  categories:  the 
first  category,  which  we  call  totally  nonideal  explosives  comprises  explosives  in  which  the  fuel 
and  oxidizer  each  appear  only  in  separate  grains.  Examples  are  AP/wax  and  AP/Al,  without  any 
energetic  binder.  It  must  be  emphasized  that  such  mixtures  are  not  practical  for  service  use,  but 
they  have  been  tried  out  as  experimental  mixtures  and  their  detonation  velocities  and  sensitivities 
have  been  measured  by  D.  Price,  et  al.  (The  AP/Al  mixture  fell  in  Price's  Category  2 
classification.)  Their  utility  is  that  they  are  excellent  examples  of  mixtures  for  whieh  the  gross 
heterogeneity,  diffusion  rates,  mixing  distances,  etc.  can  be  modeled  imencumbered  by  any  rapid 
reaction  of  an  explosive  molecule  containing  both  fuel  and  oxidizer.  It  is  this  category  of 
nonideal  explosives  that  has  been  emphasized  in  this  Phase  I  theoretical  effort,  which  is 
described  below  in  Sections  2  and  3. 

The  second  category,  which  we  call  somewhat  ambiguously  nonideal  explosives,  are  a 
mixture  of  an  ideal  explosive  and  at  least  one  separate  fuel  or  oxidizer,  e.g.,  PBXN-102 
(HMX/Al/binder)  and  PBXW-1 15(Q)  (RDX/AP/Al/binder).  Thus,  this  category  can  be  viewed 
as  an  admixture  of  ideal  and  totally  nonideal  explosives.  Most  practical  explosives  fall  within 
this  category,  particularly  those  used  by  the  Navy.  Model  development  for  this  category  of 
nonideal  explosive  is  left  for  a  future  Phase  II  effort. 

The  field  of  nonideal  explosives  has  been  extensively  investigated  for  many  years, 
particularly  by  the  underwater  explosives  community.  As  early  as  1918,  German  studies  showed 
that  the  addition  of  15-20%  aluminum  increases  underwater  performance.  There  exists  a  large 
number  of  Navy-developed  plastic  bonded  explosives  containing  ammonium  perchlorate  and/or 
aluminum:  PBXC-117(Q),  PBXN-102,  PBXN-103,  PBXN-105,  PBXN-109,  PBXW-107(Q), 
PBXW-1 14(Q),  PBXW-1 15  (Q),  all  of  which  are  cast;  PBXN-4,  PBXN-7  (Type  I),  all  of  which 
are  pressed;  PBXN-201,  which  is  extrudable  and  injection  moldable. 

E.  Anderson  provides  an  excellent  review  of  explosives,  and  in  particular  the  effects  of 
explosive  composition  ingredients,  charge  diameter,  failure  diameter,  and  confinement  on 
performance  of  nonideal  explosives.  The  Detonation  Symposia  provide  an  extensive  literature 
on  nonideal  detonations,  including  some  of  the  important  results  from  the  Lawrence  Livermore 
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National  Laboratory. 

Recently,  some  very  interesting  small-scale  laboratory  tests  by  P.  Miller  and  R.  Guirguis 
were  performed  comparing  product  gas  expansion  rates  against  a  column  of  water  for  ideal  and 
nonideal  aluminized  explosives.  JWL  equation  of  state  incorporating  a  time-dependent  late 
energy  release  was  adjusted  to  reproduce  the  experimental  nonideal  explosive  results.  These 
authors  also  calculated  the  effects  of  late  chemical  reactions  on  the  bubble  oscillation  and 
pressure  pulse  induced  in  water  by  nonideal  explosives.  Their  energy  release  was  governed  by  the 
rate  law  dA/dt=a(l-A,)’^^p'^®,  where  X,  p,  t,  and  a  are,  respectively,  the  progress  variable,  pressure, 
time,  and  a  constant. 

1.4.  Current  Modeling  Effort 

This  report  describes  a  preliminary  model  of  kinetic  and  thermo-hydrodynamic  processes 
that  define  the  nonideal  detonation  of  a  totally  nonideal  explosive:  an  AP/Al  mixture. 

Section  2  presents  a  thermo-hydrodynamic  description  of  reactive  explosive  compositions 
containing  aluminum  and,  in  particular,  a  description  of  the  nonideal  detonation  process  in 
ammonium  perchlorate/aluminum  and  CHNO  explosive/Al  compositions.  The  specific  energy 
relationships  for  these  reacting  compositions  are  presented  in  which  Mie-Gruneisen  equations  of 
state  are  derived,  including  the  temperature,  and  used  for  the  solid  components.  Also  the 
equations  of  steady  state  detonation  in  AP/Al  compositions  are  presented. 

In  Section  3  we  present  a  theory  describing  the  kinetic  processes  by  which  AI/AI2O3/AP 
products  react.  Key  here  is  the  application  of  the  phenomenon  of  anomalously  high  ionic 
conduction,  which  are  exhibited  by  the  P-aluminas,  through  “tunnels”  through  which  cations  can 
move  freely.  This  provides  a  plausible  explanation  of  events  in  Al/explosive  reactions. 

Statistical  mechanical  concepts  are  employed  to  provide  reaction  rate  equations  for  A1  and  AP. 

The  theoretical  results  of  Sections  2  and  3  provides  a  preliminary  physical  model  based 
on  first  principles,  which  will  be  integrated  into  a  compressible  flow  code  in  the  next  phase  of 
this  work,  which  is  to  improve,  extend,  and  validate  the  model. 

1.5.  Acknowledgement 

This  report  is  based  upon  work  supported  by  the  Office  of  Naval  Research  under  Small 
Business  Innovative  Research  (SBIR)  Program  Phase  I  contract  no.  N00014-95-C-0154. 

We  wish  to  thank  Drs.  Phil  Miller  and  Raafat  Guirguis  of  the  Naval  Surface  Warfare 
Center  Indian  Head  Division,  White  Oak  for  providing  profitable  discussions  and  information  on 
their  experimental  and  theoretical  work  on  nonideal  detonations  of  aluminized  explosives.  Dr. 
Jagadish  Sharma,  also  of  the  NSWCIHD,  White  Oak,  for  fmitful  discussions  on  the 
crystallography  of  ammonium  perchlorate,  Erwin  Anderson  of  NSWCIHD,  White  Oak,  for 
informative  discussions  on  nonideal  explosives,  and  to  Dr.  William  Davis,  formerly  of  the  Los 
Alamos  National  Laboratory,  for  steering  us  to  the  work  of  E.  O’Connor. 


4 


ENIG  TR  95-1 


2.  THERMO-HYDRODYNAMICS  OF  REACTING  EXPLOSIVE  MIXTURES 
CONTAINING  ALUMINUM  (Al) 

In  this  section  of  the  report,  we  formulate  a  global  thermodynamic  description  of  reacting 
explosive  compositions  containing  Al.  Such  a  description  is  required  in  order  to  perform 
analytical  and  numerical  studies  of  the  role  played  by  Al  in  the  nonideal  detonation  process  in 
ammonium  perchlorate  (AP)/A1  compositions  as  well  as  CHNO  explosive/ Al  compositions. 

In  Sections  2.1  and  2.2,  we  construct  the  specific  energy  e  relationships  for  reacting 
AP/Al  and  CHNO/Al  compositions  with  their  solid  components  governed  by  generalized  Mie- 
Gruneisen  (MG)  equations  of  state.  These  constructions  must  be  based  on  the  mass  balance 
equations  for  the  decomposing  explosive  and  the  Al  reactions  because  the  Al  combines  with 
oxygen-containing  products  of  the  explosive  to  form  aluminum  oxide  AI2O3.  We  also  derive  the 
mass  balance  equations  for  Al  spheres  reacting  with  oxygen  to  form  an  outer  coating  of  AljOj. 

In  Section  2.3,  we  presented  the  particular  MG  EOS  that  were  constructed  for  AP,  Al,  and 
AI2O3,  and  the  EOS  used  for  the  explosive’s  products. 

In  Section  2.4,  we  consider  steady-state  (SS)  detonation  waves,  and  the  relationship 
between  the  SS  reaction  zone  length  and  reaction  time  in  AP.  We  also  develop  a  protot5q)e 
model  for  SS,  nonideal  detonation  in  a  CHNO/Al  composition  to  provide  a  better  understanding 
of  the  nonideal  detonation  process  in  explosive  compositions  containing  Al,  and  demonstrate 
why  global  kinetic  expressions  must  be  constructed  for  the  Al  reactions.  A  procedure  for 
constructing  such  a  global  kinetic  expression  will  be  presented  in  the  next  Section  3  of  this 
report. 


2.1.  A  Specific  Energy  Relationship  for  a  Reacting  AP/Al  Mixture 

Because  of  the  complex  nature  of  reactive  explosive  mixtures  in  the  10-400  kbar  region 
treatments  of  detonation  phenomena  must  be  based  on  assumptions  regarding  the  significant 
chemical  and  physical  processes.  Here,  our  formulation  of  an  e-relationship  for  AP/Al  mixtures 
with  compressible,  condensed  components  is  based  on  the  following  assumptions: 

A.l .  The  AP  decomposes  to  form  a  mixture  of  products  with  a  fixed  composition 
containing  oxygen  (O2),  nitrogen  (N2)5  water  (H2O),  and  chlorine  (CI2). 

A.2.  The  AP  produces  enough  O2  to  react  with  all  the  Al  in  the  AP/Al  mixture. 

A.3.  The  AI2O3  produced  by  the  reaction  of  Al  and  O2  can  be  either  solid  a-Al203  or  y- 
AI2O3 ,  depending  on  the  thermodynamic  state. 

A.4.  Each  constituent  in  the  reacting  AP/Al  mixture  has  the  same  pressure  p  and  is 

governed  by  its  own  equation  of  state  (EOS)  relating  e,  p,  and  specific  volume  v. 


It  follows  from  A.1  and  A.3  that  our  relation  for  e  must  contain  two  reaction  coordinates, 
say  Al  and  A2,  to  account,  respectively,  for  the  decomposition  of  the  AP  and  the  reaction  of  Al 
and  O2  to  form  solid  AI2O3.  In  addition,  we  must  account  for  the  formation  of  either  a-Al203  or 
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y-AljOj  and  formulate  a  condition  to  determine  when  one  is  formed  and  not  the  other.  Such  a 
condition  will  be  formulated  later  when  we  consider  the  EOS  of  AI2O3.  At  this  stage,  it  should 
also  be  noted  that  we  will  not  consider  melting  of  either  the  A1  or  the  solid  AljOj. 

When  mixture  rules  for  the  total  internal  energy  E^  and  volume  of  such  an  AP/Al 
composition  are  specified,  these  assumptions  allow  us  to  formulate  an  e-relationship  for  this 
composition  in  terms  of  p,  the  specific  volumes  of  its  components,  and  the  reaction  coordinates 
A,,  and  Aj.  Our  formulation  of  such  an  e-relationship,  here,  will  be  carried  out  in  two  steps.  In 
the  first  step,  we  will  define  the  reaction  coordinates  Aj  and  Aj  and  derive  equations  for  the 
masses  of  the  components  in  our  reacting  AP/Al  mixture  in  terms  of  A,  and  A2.  In  the  second 
step,  we  will  use  the  mixture  rules  for  E„  and  V„,  and  the  e=e(p,v)  EOS  of  the  components  to 
derive  our  equation  for  e. 

2.1.1.  Mass  Balance  Equations  for  an  AP/Al  Composition 

We  consider  an  AP/Al  composition  with  an  initial  mass  Mq,  an  initial  mass  fraction  of 
AP,  =  m^  /Mq,  and  an  initial  mass  fraction  of  Al,  a^i  =  m^i  /Mq,  and  also  recall  that  it  is 
customary  for  such  a  mixture  to  set  Mo=100g.  To  write  the  mass  balance  equations  for  such  a 
mixture  subject  to  A.l  and  A.3,  we  must  consider  the  following  chemical  reactions: 

NH4C104  =  02  +  2H20  +  ‘/2N2  +  l/2Cl2  (R.l) 

2A1  +  3/2  O2  =  (AI2O3)*,  i=l  or  2  (R.2) 


with  i=l  and  i=2  used,  respectively,  to  denote  a-Al203  and  Y-AI2O3.  We  first  let  Mj  denote  the 
molecular  weight  of  a  species  Sj,  and  see  from  (R.1)  and  R.2)  that  m^  /M^p  ^  3mAi  /4Mai  is  the 
condition  for  our  composition  to  satisfy  A.2.  We  next  let  mj  denote  the  mass  of  a  species  S;  and 
write  the  law  of  definite  proportions  for  (R.l)  and  (R.2)  with  extents  of  reaction  and  ^2  as 


dm 


AP 


M 


AP 


dm^ 


dm. 


dm 


CL 


V2M. 


CL 


dm 


N. 


ViM, 


N, 


dm 


Al 


2M 


Al 


dm, 


dmli 


ALO^ 


=  dL 


with 


(2.1.1) 


(2.1.2) 


dm^,  =  dmn  +  dml 


(2.1.3) 


Here  the  superscripts  1  and  2  refer  to  the  respective  equations,  R.1  and  R.2. 


Formally  integrating  these  equations  and  introducing  the  reaction  coordinates 

Af.p  2M., 

A,  =  and  A,  = 


m 


AP 


m 


Al 
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and  the  mass  fractions 


2M„ 


M. 


=  a 


1  ’ 


M. 


=  a. 


2M. 


M. 


=  a. 


Cl. 


2M,,  3Af 


2M. 


=  a. 


M. 


=  Pi  . 


2  _ 


2M. 


ylf 

gives  the  equations  for  the  masses  of  the  species  in  our  reactive  AP/Al  mixture  as 

~  ^Ai  P2VP1 

^Ai  =  -^2) 

^ALO,  ~  ^^Al  /Pl)('^2~^3) 


m 


H.O 


m 


0  ^ 


m.T  ~  nt 


0  T 

A 


AP'^3  I 


^Cl,  =  ^ap%K 


(2.1.4) 

(2.1.5) 

(2.1.6) 

(2.1.7) 


=  P2 


(2.1.8) 

(2.1.9) 

(2.1.10) 


with  m^i  Q  used  to  denote  the  mass  of  (Al203)s.  It  follows  from  these  equations  that  the  masses 
of  the  condensed  components  are  defined  by  Eqs.  (2.1.4),  (2.1.9),  and  (2.1.10)  and  that  the  mass 
of  the  fluid  component,  mf,  is  obtained  by  summing  Eqs.  (2.1.5)  through  (2. 1 .8).  We  will  now 
present  the  mixture  rules  and  the  EOS  used  in  the  formulation  of  an  e-relationship  for  our 
reacting  AP/Al  mixture. 


2.1.2.  Mixture  Rules  and  Equations  of  State 

Here,  because  it  is  convenient  to  express  mixture  rules  in  terms  of  extensive  variables,  we 
write  En,  and  V„  for  our  AP/Al  composition  in  terms  of  the  extensive  internal  energies  and 
volumes  of  its  components  as 

'■=!  or  2  (2.1.12) 

K  --  '=1  O'-  2  (2.1.13) 

In  this  case,  the  specific  internal  energy  and  volume  of  the  mixture  are  related  to  En,  and  by 
the  equations,  e=En,/Mo  and  \=VJMq.  We  now  introduce  the  EOS  used  to  calculate  the  energies 
on  the  right-hand  side  of  Eq.  (2.1.12). 

We  write  the  e=e(p,v)  EOS  for  the  condensed  components  in  Gruneisen  form  as 

^AP  =  ^  +  ^t^’sAP^y)  (2.1.14) 

^  AP 


P^Al 


^  <SAfy) 


(2.1.15) 


(A/,;)' 


P'^'i  i 

'f’Alfi^  ^  ^  Vog.(E) 


/=1  or  2 


(2.1.16) 
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In  writing  these  equations,  we  have  used  (Ah®);  to  denote  the  specific  heat  of  formation  of 
a  species  Sj  and  have  reduced  their  number  by  using  i  as  a  subscript  or  a  superscript  with  the 
values  1  and  2,  respectively,  to  denote  a-Al203  and  y-AljOj.  In  addition  the  g  terms  are  written 
as  functions  of  the  reduced  volume  V=v/vo,  with  the  understanding  that  V  in  gi(V)  denotes 
V=(v/vo)i.  The  energy  terms  for  the  condensed  components  in  Eq.  (2. 1 .12)  can  readily  be 
obtained  by  multiplying  Eqs.  (2.1.14),  (2.1.15)  and  (2.1.16)  by  the  expressions  for  their 
respective  masses  given  by  Eqs.  (2. 1 .4),  (2. 1 .9),  and  (2.1.10). 

We  write  the  Ef=E(p,Vf,  mf)  EOS  for  the  fluid  as 

~  ^  ~  ’  (2.1.17) 

where  the  reduced  volume  V  in  gf  at  this  stage  need  not  be  defined.  Only  the  specific  heat  of 
formation  for  H2O  need  be  included  in  Eq.  (2.1.17)  because  the  heats  of  formation  of  all  the  other 
fluid  constituents  are  zero.  We  can  now  derive  a  specific  energy  relationship  for  our  AP/Al 
mixture  as  follows. 


2.1.3.  Derivation  of  the  Specific  Energy  Relationship 


We  first  derive  an  expression  relating  En,  and  by  substituting  the  equation  for  Ef 
obtained  by  eliminating  Vf  between  Eqs.  (2.1.13)  and  (2.1.17)  into  Eq.  (2.1.12).  The  result  of 
these  operations  leads  to  the  equation 


•  (2.1.18) 


We  then  multiply  Eqs.  (2.1.14)-(2.1.16)  for  the  specific  energies  e^p,  e^i,  and  e^,  ^  of  the 
condensed  components  by  the  expressions  for  their  respective  masses  in  Eqs.  (2.1.4).  (2.1.9)  and 
(2.1.10)  to  obtain  expressions  for  E^p,  E^i,  and  E^,  ^  ,  substitute  these  expressions  into  Eq. 
(2.1.18),  and  collect  similar  terms.  The  equations  obtained  in  this  way  for  the  heat  of  formation 
terms  denoted  by  E^,  the  pressure  terms  denoted  by  Ep,  and  the  remaining  volume  terms  denoted 
by  Eg  are  given  below. 


Collecting  the  specific  heat  of  formation  terms  for  the  components  leads  to  the  following 
equation: 


Pi 


(2.1.19) 


We  now  let  E°  and  Q,  denote,  respectively,  the  heat  of  formation  of  our  AP/Al  composition  and 
the  chemical  energy  liberated  when  1  g  of  AP  decomposes  according  to  R.  1 .  We  also  let  Qj 
denote  the  chemical  energy  liberated  when  1  g  of  A1  reacts  to  form  a-Al203  (i=l)  or  Y-AI2O3  (i=2) 
according  to  R.2.  We  then  write  Eq.  (2.1.19)  in  terms  of  these  parameters  as 

.  (2,1.20) 

Collecting  the  pressure  terms  leads  to  the  equation 
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^nt  A  Ap  f.  AJ  n  A 1 

Ep  =  p—^ — - — c^-KKp-^p—^^Aii^-KKi  + 

where  =  (n/T^p-\),  =  (n/r^,-l),  and  A  =  (n/F-l)  . 

Collecting  the  remaining  volume  terms  leads  to  the  equation 


(2.1.21) 


+  m^,(l-X^)[VogJV)^v^^^(F)] 
+  /PMVog.(^+V^^  (V)] 


(2.1.22) 


The  equation  for  e  for  our  mixture  is  obtained  by  adding  Eqs.  (2.1.20),  (2.1.21),  and 
(2.1.22)  and  dividing  by  Mg.  Because  of  space  limitations,  it  serves  our  purpose  here  to  set  gf=0, 
and  n=k-l  and  only  present  the  energy  relationship  when  the  products  are  assumed  to  be 
polytropic.  In  this  case,  the  e-relationship  can  be  written  as 

(2.1.23) 

■^^ApC^  -^)gAp(^^^Al(l  -^2)^^,(^+«>2g,(T)  , 
where  -  cc^pQj,  ?2  “  ^apQz’  -^ap  ~  ^ap^ap’  -^ai  ~  ^a/^ai’-^i  ~  ^At^/Pi> 

0  „0  AP  0  0  Al  J  0  0  1,0 

^AP  ~  ^AP^O  »  “  ^Al^O  5  ^d  a,.  -  a^^Vg/Pj. 


At  this  stage,  Eq.  (2. 1 .23)  cannot  be  used  to  calculate  the  specific  energy  of  our  reacting 
AP/Al  mixture  because  we  have  no  prescription  for  calculating  the  specific  volumes  of  the 
condensed  components  when  p,  v.  A,,,  and  Aj  are  known.  A  prescription  to  address  this  difficulty 
will  be  presented  after  we  have  derived  the  mass  balance  equations  for  reacting  Al  particles  in  the 
next  section. 


2.1.4.  The  Mass  Balance  Equations  for  Al  Spheres  Reacting  with  Oj  to  Form  an 
Outer  Coating  of  AljOj  and  the  Reaction  Coordinate  A2. 


We  consider  a  shocked  sphere  of  Al  with  an  initial  radius  rjj,  specific  volume  v”,  and 
mass  w?.  The  initial  mass,  specific  volume  and  radius  of  this  Al  sphere  are  then  related  by  the 
equation 


=  —  •  (2.1.24) 

3vj 

We  then  assume  that  a  mass  fraction  of  Al,  Aw°,  has  reacted  with  the  O2  in  the 
surrovmdings  to  produce  a  sphere  of  radius  R.  This  sphere  consists  of  an  outer  shell  of  AljOj 
with  mass  Wj,  specific  volume  Vj  and  a  thickness  /*,  that  covers  an  inner  core  of  Al  with  a  radius 
R-A  In  this  case,  the  mass  of  Al  in  the  inner  core,  w,,  satisfies  the  equations 
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4%(R-l‘f  4lTr^ 
w.  -  - ^  =  - -(1-A)  , 


1 


1 


3v 


H 

1 


which  gives  the  following  equation  for  X: 


1 


H 


(2.1.25) 


(2.1.26) 


We  note  that  Eq.  (2.1.26)  satisfies  the  conditions,  A,h=0,  when  /-O  and  R=rH,  and  A=1  when  the 
thickness  of  the  AI2O3  coating  equals  the  radius  of  the  sphere.  We  now  use  the  mass  balance 
equation  from  R.2  to  derive  another  equation  for  R  and  thereby  obtain  the  relationship  between  t 
and  X.  We  accordingly  integrate  the  equation  obtained  by  setting  mA,=W],  Ma,=Mi,  w  ^  =wj, 
and  M2  in  Eq.  (2.1.2)  to  obtain  the  following  equations  relating  w,  and  W2: 


w?-Wi  =  2Mi^2  (2.1.27) 

w1  =  M2^2  (2.1.28) 


Eliminating  ^2  between  Eqs  (2.1.27)  and  (2.1.28)  then  gives  the  equation 


Awq  =  (2M1/M2)  W2, 


(2.1.29) 


which  we  rewrite  as 


■j  ^  M  3  ' 

3vj  •'“2  ^  Vj 


(2.1.30) 


The  combination  of  Eqs.  (2.1.26)  and  (2.1.30)  then  leads  to  the  following  equation  relating  R  and 
X: 


TT  ^  TT  ^ 


1/3 


2r, 


(2.1.31) 


where  V,=ViMi  denotes  the  molar  volume  of  A1  and  V2=V2,  M2  denotes  the  molar  volume  of 
AI2O3.  Setting  A=1  in  Eq.  (2. 1.31)  gives  the  final  radius  of  the  sphere  when  all  the  A1  has  reacted 
to  form  AI2O3.  The  combination  of  Eqs.  (2.1.26)  and  (2.1.31)  then  gives  the  equation  relating  /' 
and  X  as 


/ '  vl  V 

—  -  [—X  +  — L(1-A)]‘^^ 
2vf  v" 


(-^)'''(1-^) 
Vf 


1/3 


(2.1.32) 


Having  introduced  the  reaction  coordinate  X  for  a  single  A1  particle,  we  must  now  relate 
this  coordinate  to  the  global  coordinate  A2  defined  by  the  equation  mAi=mAi(l-A2).  For  simplicity, 
we  will  assume  that  our  AP/Al  mixture  contains  N^i  particles  of  A1  with  the  same  initial  radius  ro, 
and  write  the  equation  for  m^i  as  mAi=NA,4TT:ro/3v^'.  We  now  multiply  Eq.  (2.1.25)  by  N^i  to 


10 


ENIGTR95-1 


obtain  the  total  amount  of  unreacted  A1  in  our  mixture,  when  the  covering  layers  of  AI2O3  have  a 
thickness  T,  as 

AlzN.rl 

^Ai  =  ^Ai^i  ^  ^  (1"^)  •  (2.1.33) 

3v" 

The  mass  balance  condition  (j-Qlvf)  =  (rpvf)  then  allows  us  to  rewrite  Eq.  (2.1.33)  as 
mA|=mAi(l-A)  and  show  that  X=X2  when  all  the  A1  particles  in  our  AP/Al  mixture  are  assumed  to 
have  the  same  size. 

2.1.5.  Implementation  of  the  e-Relationship. 

Our  prescription  to  make  Eq.  (2.1.23)  into  a  working  relationship  is  based  on  the  premise 
that  the  components  of  our  reacting  AP/Al  mixture  do  not  attain  thermal  equilibrium  during  the 
detonation  process.  Consequently,  because  an  explicit  treatment  of  heat  conduction  is  beyond 
the  scope  of  Phase  I  tasks,  addtional  assumptions  must  be  made  about  the  transfer  of  heat  among 
these  components.  Here,  we  will  treat  the  AP  and  A1  in  a  similar  way,  and  assume  that  no 
appreciable  amount  of  heat  is  transferred  into  either  the  interior  of  a  burning  AP  grain  or  the 
interior  of  a  reacting  A1  particle,  as  the  reacting  AP/Al  mixture  expands.  In  this  case,  a  grain  of 
AP  and  a  particle  of  A1  entering  the  shock  are  thereafter  constrained  to  their  respective  isentropes 
determined  by  their  thermodynamic  states  attained  at  the  shock  firont.  In  contrast  to  this 
treatment  of  AP  and  Al,  we  assume  that  the  heat  transfer  process  between  AljOj  and  the  fluid 
component  as  so  efficient  that  the  AI2O3  coatings  on  the  Al  particle  attain  thermal  equilibrium 
with  the  decomposition  products  of  the  AP.  The  consequences  of  these  assumptions,  which 
make  Eq.  (2. 1 .23)  a  workable  e-relationship  will  be  considered  later. 

2.2  A  Specific  Energy  Relationship  for  a  Reacting  CgH^NcOyAl  Mixture. 

Our  formulation  of  an  e-relationship  for  a  CaH^NpO^Alf  composition  follows  a  procedure 
similar  to  that  used  in  our  formulation  of  the  e-relationship  for  AP/Al  mixture.  To  be  more 
specific,  we  adopt  the  assumption  A.3,  apply  assumption  A.4  to  the  CaH^NgOdAlf  system,  and 
make  the  following  assumptions  about  the  decomposition  of  the  explosive  CaHbN(,Od,  and  the 
reactions  of  its  products  with  the  Al: 

A.  5  The  CaH^NjOd  explosive  decomposes  to  give  products  with  a  fixed  composition,  as 
prescribed,  by  Kamlet*. 

A.6  The  Al  reacts  with  the  CO2  and  H2O  produced  by  CaMbN^Oj,  and  the  H2  liberated 
from  the  H2O  reacts  with  the  C  in  the  explosive’s  products  to  form  hydrocarbons. 

2.2.1.  Mass  Balance  Equations  for  a  CaH|,NaOjAlf  Composition 

We  again  consider  a  mass  Mo=100g  of  explosive  mixture  with  an  initial  mass  fraction  of 
X,  a°=m^/Mo,  and  an  initial  mass  fraction  of  Al,  aAi“rnAi/Mo.  It  is  clear  for  this  X/Al  mixture  that 
m®|=fMA|  and  that  m°=Mo-fMAi.  Following  Kamlet,  we  write  the  equation  for  the  decomposition 
of  our  explosive  component  as 
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CaHbNjOjAlf  bjH20+aiC02+C|N2+a2C"i'fAl 


(R.3) 


with  bi=b/2,  ai=d/2-b/4,  Ci=c/2,  and  a2=a-(d/2-b/4).  To  consider  the  heat  produced  by  (R.3),  we 
denote  the  heat  of  formation  of  a  mole  of  a  species  S  as  (AH°f)s,  and  write  the  standard  heat  of 
decomposition  of  X  in  the  composition,  Q^,  as 


(2.2.1) 


and  its  specific  heat  of  decomposition  as 


X 


-M(LhX] 


(2.2.2) 


because  Mx=m°. 

We  now  consider  the  reactions  of  A1  with  the  explosive’s  products.  We  write  the 
equation  for  the  A1  in  the  composition  reacting  with  the  CO2  as 

fAl*^CO^  .  .  (R.4) 

Examination  of  (R.3)  and  (R.4)  then  shows  that  all  the  A1  in  the  composition  can  react  with  the 
CO2  when  ai  >3E2.  We  write  the  equation  for  the  A1  in  the  composition  reacting  with  the  H2O  as 
and  similarly  see  that  the  condition  for  all  the  A1  to  react  with  the  H2O  is  bi^3f/2.  Because  (R.5) 


produces  Elj,  we  must  also  consider  the  reaction  of  the  H2  with  the  C  produced  by  the  explosive. 
For  simplicity,  we  only  consider  one  reaction  between  H2  and  C  and  write  this  reaction  as 


(R.6) 


Examination  of  (R.3)  and  (R.6)  then  shows  that  the  condition  for  all  the  hydrogen  produced  by 
the  H2O  to  react  with  the  C  is  a2>f.  On  the  added  assumption  that  the  E12  produced  by  the  H2O 
reacts  immediately  with  the  C  produced  by  the  explosive,  we  combine  (R.5)  and  (R.6)  and  write 
the  equation  for  this  AI-H2O-C  reaction  as 


and  the  equation  for  its  standard  heat  of  reaction,  Q^i,  as 

q1,  =  ] 


(R.7) 


(2.2.3) 


To  be  more  specific  about  such  X-Al  compositions,  we  consider  a  composition  similar  to 
H6  with  a=1.8869,  b=2.6123,  c=1.5796,  d=1.9667,  and  f=0.7827.  Because  ai=0.3303  and 
3fi2=1.1740,  we  see  that  3f/2>a,  for  this  composition,  and  X  does  not  produce  enough  CO2  to 
consume  all  the  Al.  Alternatively,  because  bi=l  .3061>3f/2,  and  a2=l  .5566>f,  the  H2O  produced 
by  X  can  consume  all  the  Al,  and  the  C  produced  by  X  can  react  with  all  the  H2  product  from  the 
H2O  by  (R.3).  For  these  reasons,  it  is  convenient  here  to  neglect  the  reaction  of  Al  and  CO2  in 
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this  composition  and  only  consider  the  two  global  exothermic  reactions,  (R.3)  and  (R.7). 

In  order  to  account  for  (R.3)  and  (R.7)  in  the  detonation  process,  we  introduce  the  extent 
^3  for  (R.3),  the  extent  of  reaction  for  (R.7),  and  write  the  law  of  definite  proportions  for  these 
reactions  as 


dm  dm 


H,0 


M  bM. 


dm, 

cM 


dm 


CO. 


^\^co. 


dm. 


(2.2.4) 


dm, 


dm 


Al 


dm 


H.O 


dm , 


dm. 


fMr 


Al 


(2.2.5) 


with 


dm. 


dmt  n  +  dml 


(2.2.6) 


dm„  =  dm^  +  dml 


(2.2.7) 


Formally  integrating  these  equations,  setting  ^3=A,3  and  ^7=^7,  and  introducing  the  mass 
fractions 


c,M 


«i  = 


r'-'AT, 


m. 


bM. 


,  ct,  = 


,  a. 


^i^co. 


2  „  _ 


.  a.  = 


m. 


m. 


m. 


allows  us  to  write  the  equations  for  the  masses  of  the  species  in  our  reactive  X/Al  mixture  as 


m 

m 

m 

m 


(2.2.8) 

N,  =  ““1^3 

(2.2.9) 

cOj  “  m^a^X^ 

(2.2.10) 

(2.2.11) 

c 

m^a^X^  -  fM^X,j 

(2.2.12) 

=  —M  X 

(2.2.13) 

Al 

-^7) 

(2.2.14) 

'ai,o. 

(2.2.15) 

The  Eqs.  (2.2.8),  (2.2.14),  and  (2.2.15),  respectively,  give  the  masses  of  the  condensed 
components  X,  Al,  and  AI2O3  in  terms  of  the  reaction  coordinates,  and  the  equation  for  the  mass 
of  the  other  species  can  be  obtained  by  summing  Eqs.  (2.2.9)  through  (2.2.13)  as 


(2.2.16) 


To  simplify  our  treatment  here,  but  accoimt  for  the  presence  of  solid  C  in  the  detonation 
products,  we  make  the  addition2il  assumption  that  the  species  CO2,  H2O,  N2,  and  C2H6  form  a 
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polytropic  fluid  mixture  with  an  index  k  and  that  C  is  a  polytropic  solid  with  the  same  index.  We 
can  then  proceed  with  our  formulation  of  an  e-relationship  for  our  X/Al  mixture. 

2.2.2.  Mixture  Rules  and  Equations  of  State  (EOS) 

We  again  use  and  to  denote  the  internal  energy  and  volume  of  our  X/Al 

composition,  and  express  E^  and  in  terms  of  the  internal  energies  and  volumes  of  its 
components  as 

^n,  =  (2.2.17) 

We  then  write  the  e=e(p,v)  equation  of  state  (EOS)  of  the  explosive  as 


X 


(2.2.18) 

(2.2.19) 


use  Eqs.  (2.1.15)  and  (2.1.16)  to  describe  A1  and  AI2O3,  and  write  the  Ef=E(p,Vf,  mj)  EOS  for  our 
products’  component  as 


EIl 

Y-l 


(2.2.20) 


with  j=C02,  H2O,  C2H6,  and  ,  and  ^  defined  by  Eqs.  (2.2.10),  (2.2.1 1),  and 

(2.2. 1 3).  Our  derivation  of  a  specific  energy  relationskip  for  the  X/Al  mixture  now  parallels  that 
for  the  AP/Al  mixture. 


2.2.3.  Derivation  of  the  Specific  Energy  Relationship 


We  first  derive  an  expression  relating  E„  and  by  substituting  the  equation  for  Ef, 
obtained  by  eliminating  Vf  between  Eqs.  (2.2.18)  and  (2.2.10),  into  Eq.  (2.2.17).  These 
operations  give 


pVa, 

E_  =  (£ -^)  +  (£,r^)  +  (<o. 


Y-l 


Y-l 


P^ko.  „  *  0 

- - ).  +  Hm  {Lhf).+ - 

Y-l  '  '  Y-l 


(2.2.21) 


We  then  multiply  Eqs.  (2.2. 1 9),  (2. 1 . 1 5),  and  (2.1.1 6),  for  the  specific  internal  energies  of  the 
condensed  components,  by  the  expressions  for  their  respective  masses  given  in  Eqs.  (2.2.8), 
(2.2.14),  and  (2.2.15),  to  obtain  expressions  for  E,;,  E^i,  and  q  ).,  substitute  these  expressions 
into  Eq.  (2.2.21),  and  collect  similar  terms.  The  resulting  equations  for  the  heat  of  formation 
terms  denoted  by  E^,  the  pressure  terms  denoted  by  Ep,  and  the  remaining  volume  terms  denoted 
by  Eg  are  given  below. 

The  equation  for  E^.  is  obtained  as 

=  <-Q.K-Qlik  ’  (2.2.22) 

with  E°=m°(Ahf),j  because  (Ahf)Ai=0.  The  equation  for  Ep  is  obtained  as 
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/, 


(2.2.23) 


with  A^- ^  ^  -1  ,  ^  -1  ,  a«c?  defined  previously  for  the  AP/Al 

^  X  ^  Al 


r 


composition.  The  equation  for  E  is  obtained  as 


E,  =  r;(i -x,)g^(V)*Ffo -x,)gjn*^M„oyjx,g^n  ,  (2,2.24) 

with  r/=w^“v/  and  . 

Finally,  adding  Eqs.  (2.2.22),  (2.2.23),  and  (2.2.24)  together  and  dividing  the  resulting 
equation  by  Mq,  gives  the  e-relationship  for  our  X/Al  mixture  as 


with 


e-e  = 


^ 

+  [«;(!  -X^)g,(V)+a^i(l  -k^)gJV)  V.Wl 


(2.2.25) 


93  =  .  ^ArSA/K  »  ^x  =^,^x  . 

Having  derived  e-relatinships  for  AP/Al  and  X/Al  mixtures  in  terms  of  general  MG  EOS,  we 
must  now  formulate  the  particular  EOS  that  are  needed  for  their  implementation. 


2.3  Equations  of  State  for  the  Components  in  AP/Al  and  X/Al  Compositions. 

Because  A1  is  present  in  both  of  our  explosive  compositions,  we  will  formulate  an  EOS 
for  A1  before  considering  the  EOS  of  the  other  components. 


2.3.1.  A  Complete  Equation  of  State  for  A1 

Difficulties  arise  in  modeling  the  reactions  of  a  metal  with  an  explosive’s  products  behind 
the  detonation  front  when  the  EOS  of  the  metal  is  based  on  shock  wave  data  and  is, 
consequently,  valid  only  for  values  of  the  specific  volume  v,  v<Vo=Po',  where  p  denotes  the 
density  and  the  subscript  0  denotes  the  unshocked  state.  An  approach  to  resolve  this  difficulty, 
which  is  different  from  the  one  presented  here,  has  been  formulated  by  J.W.  Enig^.  To  be  more 
specific  about  this  difficulty,  we  consider  an  e=e(p,v)  EOS  based  on  the  assumption  that  the 
Gruneisen  parameter  F  is  only  a  function  of  v.  (The  thermodynamic  conditions  for  this 
assumption  to  be  valid  will  be  presented  later.)  In  this  case, 

(0e/ap),  =  v/r(v),  (2.3.1) 
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and  the  e=e(p,v)  relation  follows  by  formally  integrating  Eq.  (2.3.1),  yielding 

e  =  pv/r(v)+g(v),  (2.3.2) 

where  g(v)  is  an  arbitrary  function  of  v.  When  r(v)  is  known,  the  knowledge  of  eR(v)  along  a 
reference  curve  Pr=  Pr(v)  in  the  (v-p)  plane  allows  us  to  determine  g(v)  over  the  volume  range, 
v,:sv^v2,  spanned  by  Pr(v)  and  to  rewrite  Eq.  (2.3.2)  as 

e-eR(v)  =  (p-pR(v))v/r(v)  for  v,  ^  vsvj.  (2.3 .3) 

Because  shock  wave  studies  of  a  material  provide  its  Hugoniot  curve  in  the  (v-p)  plane, 

Ph  =  Ph(v,Vo)  for  v^vo,  (2.3.4) 

and  because  the  energy  change  along  a  Hugoniot  curve  is  defined  by  the  Hugoniot  relation 

2(eH-eo)  =  Ph(vo-v),  (2.3.5) 


it  is  convenient  to  write  an  e=e(p,v)  EOS  calibrated  with  shock  wave  data  as 

e-en  ==  (p-PH)v/r(v).  (2.3.6) 

When  writing  an  e=e(p,v)  EOS  in  this  form,  one  must  remember  that  its  validity  is  limited 
strictly  to  the  volume  range,  v^Vq,  sparmed  by  the  Ph=Ph(v,Vo)  Hugoniot  curve. 

To  provide  more  background  for  our  EOS  task,  we  now  discuss  Debye’s  approach  to 
solids  to  establish  thermodynamic  conditions  for  F  to  be  a  function  of  v. 

2.3.1. 1.  Thermodynamics  of  a  Debye  Solid 

We  let  T  and  Uo(v)  denote,  respectively,  the  absolute  temperature  and  the  specific  internal 
energy  along  the  zero  degree  (T=0  K)  isotherm,  and  write  the  equation  for  the  specific  Helmholtz 
free  energy  f  as 

f  =  Uo(v)+TF(0/T),  (2.3.7) 

where  the  Debye  temperature  0  is  assumed  to  be  a  function  only  of  v  and  is  related  to  the 
maximum  frequency  that  can  be  propagated  in  the  lattice  0)^,  by  the  equation  0=(h/k)(o^/2Tr. 
Equation  (2.3.7)  is  used  to  derive  an  equation  for  F/v,  rewritten  as 

F/v  =  (ap/aT)yc„  (2.3.8) 

where  denotes  the  specific  heat  at  constant  volume. 

To  obtain  the  equation  for  (ap/5T)y,  we  first  differentiate  Eq.  (2.3.7)  partially  with  respect 
to  V  to  obtain  the  following  equation  for  p: 
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P  =  Po(v)  -F'd0/dv  (2.3.9) 

where  po(v)=-duo/dv  denotes  the  pressure  along  the  zero  degree  isotherm  and  F'=dF/d(©/T).  We 
then  differentiate  Eq.  (2.3.9)  partially  with  respect  to  T  to  obtain  the  equation  for  (dp/dT)^  as 

(dp/dX)^  =  (0F"/T^)d0/dv  (2.3.10) 

with  F"=d^F/d(0/T)^.  To  obtain  the  equation  for  Q,  we  differentiate  Eq.  (2.3.7)  twice  with 
respect  to  T  to  obtain  the  following  equations  for  the  specific  entropy  s  and  its  derivative 
{ds/dTX, 


s  =  -F+(0/T)F'  (2.3.11) 

(ds/dT),  =  -(0VT')F"  (2.3.12) 

and  we  then  use  the  identity  T(ds/dT)y=Cy  to  obtain  the  equation  for  Q  as 

C=-(0/T)2F"  (2.3.13) 

The  combination  of  Eqs.  (2.3.8),  (2.3.10),  and  (2.3.13)  then  gives  the  equation  for  E  as 

r  = -(v/0)d0/dv,  (2.3.14) 

which  shows  that  F  is  a  function  only  of  v. 


The  assumption  that  E  is  a  function  of  v  is  thus  valid  when  the  contribution  of  the  lattice 
vibrations  to  f  are  realistically  accounted  for  by  the  dependence  of  the  Debye  temperature  0  on  v. 

To  complete  our  thermodynamic  discussion  of  EOS  with  E=E(v),  we  combine  the 
identity  e=f+Ts  and  Eqs.  (2.3.7),  (2.3.9),  (2.3.1 1),  and  (2.3.14)  to  derive  the  e=e(p,v)  EOS  for  a 
Debye  solid  as  follows.  We  first  substitute  Eqs.  (2.3.7)  and  (2.3.1 1)  into  the  identity  for  e  to 
obtain  the  equation 

e=Uo(v)  +  0F',  (2.3.15) 

which  is  the  complete  e=e(s,v)  EOS  for  the  solid  because  0  is  only  a  function  of  v  and  s  is  only  a 
function  of  0/T.  We  combine  Eq.  (2.3.1 1)  with  the  Eqs.  (2.3.9)  and  (2.3.4)  to  obtain  the  e=e(p,v) 
relation  as 

e  =  Uo(v)  +  (p-po(v))v/E  (2.3.16) 

Comparing  Eqs.  (2.3.16)  and  (2.3.3),  then  shows  that  the  natural  reference  curve  in  the  (v-p) 
plane  for  the  Debye  solid  is  the  zero  degree  isotherm,  and  comparing  Eqs.  (2.3.16)  and  (2.3.2) 
shows  that  the  function  g(v)  for  the  Debye  solid  is 

g(v)  =  Uo(v)  -  0-'(dv/d0)duo/dv  (2.3.17) 
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2.3.1.2.  Approach  to  the  EOS  Task 

Our  approach  to  formulate  a  realistic  EOS  for  A1  for  detonation  calculations  was  based  on 
two  observations.  The  first  observation  is  that  to  a  good  approximation,  for  A1  subject  to  the 
ranges  of  pressures  and  temperatures  of  interest,  is  constant  and  the  Gruneisen  parameter  has  a 

constant  value  r=2.  The  second  observation  is  that  D.J.  Pastine'®  had  considerable  success 
calculating  shocked  states  in  face-centered  cubic  metals  by  using  a  Morse  potential  to  represent 
their  zero  degree  isotherms.  This  Morse  potential  can  be  written  for  a  gram  atomic  weight  of 
metal  in  terms  of  the  distance  paramter  y=l-(v/vo)''^,  with  Vq  used  to  denote  v  at  O^K,  a  parameter 
b,  and  the  latent  heat  of  sublimation  Uq  as 

U  (v/vo)  =  Uo(exp(2by)-2exp(by))  (2.3.1 8) 

Our  approach,  based  on  these  observations,  can  be  stated  as  follows.  Use  the  values  0^=0.221 
cal/g,  r=2,  and  shock  wave  data  for  Al,  and  the  Morse  potential  to  generate  reference  curves 
®r=®r(v)  and  Pr=Pr(v)  and,  thereby,  construct  an  EOS  for  Al  that  is  not  restricted  to  the  volume 
range  spanned  by  its  Hugoniot  curve  Ph=Ph(v,Vo). 

It  was  thought  necessary,  in  implementing  this  approach,  to  construct  an  EOS  for  Al  in 
the  region  spanned  by  its  Hugoniot  curve,  and  then  ascertain  the  thermodyanic  compatability  of 
this  EOS  with  the  Morse  potential.  To  provide  the  background  for  this  implementation,  we  first 
present  a  method  for  using  shock  wave  data  to  construct  a  complete  EOS  for  a  material  with 
constant  values  of  and  F.  We  then  present  our  procedure  for  testing  the  compatibility  of  this 
EOS  and  the  Morse  potential. 

After  completing  this  background,  the  application  of  our  approach  to  Al  is  given. 

2.3.1.3.  A  Complete  EOS  Based  on  Shoek  Wave  Data 

For  our  present  purpose,  it  is  convenient  to  consider  a  complete  EOS  as  an  EOS  that 
allows  us  to  calculate  the  values  of  e,  T,  and  s  at  a  point  in  the  (v-p)  plane.  In  this  case,  when 
and  r  are  known  constants,  we  claim  that  knowledge  of  the  Hugoniot  curve  Ph=Ph(v,Vo)  through 
the  initial  state  (v^,  po,  eo,  Sq,  Tq)  provides  enough  information  to  construct  a  complete  EOS  in  the 
region  of  the  (v-p)  plane  spanned  by  Ph==Ph(v,Vo)  above  the  Tg  isotherm.  We  substantiate  this 
claim  as  follows. 

For  convenience,  we  assume  that  the  Ph=Ph(v,Vo)  curve  satisfies  the  equation 

Ph  =  a2(Vo-v)/(vo-B(vo-v))2  (2.3.19) 

obtained  from  a  linear  shock  velocity  U-particle  velocity  u  relation,  U=a+Bu,  and  then  consider 
the  following  differential  equation  for  shock  temperature  Tjj, 

dTH+(r/v)THdv  =  (2Q-‘  [(vo-v)dpH+PHdv]  (2.3 .20) 

obtained  by  equating  the  differential  dey  calculated  from  the  e=e(T,v)  EOS  and  Eq.  (2.3.5). 
Formally,  integrating  Eq.  (2.3.20)  gives  the  equation  for  Tjj  as 
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*  (2C,v'')-'fv''[(v„-v)(^)tp„]iv  ,  (2.3.2 

V  J  dv 

”0 

which  shows  that  the  shock  temperature  can  be  calculated  along  the  Hugoniot  curve  p=Ph(v,Vo) 
when  Pjj  is  a  known  function  of  v,  as  exemplified  by  Eq.  (2.3.19). 

After  Tjj  has  been  calculated  as  a  function  of  v,  the  equation 


C,r(T-TH)-v(p-pH)  (2.3.22) 

can  be  used  to  calculate  temperature  along  isochores  emanating  from  the  Tq  isotherm.  Because  e 
can  also  be  calculated  along  these  isochores  from  Eq.  (2.3.6)  when  F  is  a  known  constant,  we 
need  finally,  to  derive  an  equation  for  calculating  s.  Consequently,  we  integrate  the  differential 
ds  obtained  from  the  s=s(T,v)  EOS,  to  obtain  the  following  equation  for  the  entropy  along  the 
Hugoniot  curve, 

SH-So+C,[ln(TH/To)+rin(v/vo)],  (2.3.23) 

and  write  the  equation  for  the  entropy  along  isochores  emanating  from  Ph-Ph(v,Vo)  as 

s  =  Sh  +  Cv  ln(T/TH)  (2.3.24) 

Combining  Eqs.  (2.3.23)  and  (2.3.24)  then  gives  the  equation  for  calculating  the  entropy  as 

s  =  So+C^[ln(T/To)+rin(v/Vo)]  (2.3.25) 

It  is  now  convenient  to  write  the  equations  for  the  isotherm  and  the  isentrope  through  the 
initial  state  (Vo,Po=0)  because  we  may  use  the  Morse  potential  to  try  to  match  one  of  these  curves 
later.  Setting  T=To  in  Eq.  (2.3.22)  gives  the  equation  for  the  Tq  isotherm  in  the  (v-p)  plane  as 

P  ~  Ph  ■  (Th'Tq)  CyF/v  (2.3.26) 

And  substituting  into  Eq.  (2.3.22)  the  expression  for  T  obtained  by  setting  s=S(,  in  Eq.  (2.3.25), 
gives  the  equation  for  the  isentrope  through  the  initial  state  (Vo,po=0)  as 

P  =  Ph-[Th-To(vo/v)^C^/v  (2.3.27) 

2.3.I.4.  The  Compatability  Condition  for  the  Morse  Potential  and  an  EOS 

Calibrated  with  Shock  Wave  Data 

Because  Eq.  (2.3.18)  for  the  Morse  potential  was  originally  written  for  the  OK  isotherm 
and  the  specific  volume  of  metals  at  0  K  and  p=0  are  not  really  known,  we  will  use  this  equation 
here  with  Vq  considered  to  be  the  specific  volume  in  the  unshocked  state.  In  this  case,  Uo  is  the 
latent  heat  of  sublimation  at  Tq.  We  assume  that  Eq.  (2.3.18)  represents  the  energy  along  an 
adiabat  so  that  the  equation  for  the  isentrope  passing  through  (Vq,  p=0)  is  obtained  by 
differentiating  Eq.  (2.3.18)  with  respect  to  v  as 
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p  =  -(dU  /dy)(dy/dv)  (2.3 .28) 

and  p=0  when  d  U /dy=0  and  v=Vo.  Our  compatability  condition  for  the  Morse  potential  and  an 
EOS  based  on  shock  wave  data  is  that,  for  the  pressure  range  of  interest,  we  can  find  a  value  of  b 
for  which  the  values  of  p  calculated  along  the  isentrope  defined  by  Eq.  (2.3.28)  match  those 
calculated  along  the  isentrope  defined  by  Eq.  (2.3.27). 

2.3. 1.5.  Application  of  This  Approach  to  A1 

In  order  to  determine  if  the  Morse  potential  can  be  incorporated  into  the  EOS  for  A1 
based  on  shock  wave  data,  values  of  the  pressure  calculated  for  A1  using  Eq.  (2.3.27)  must  be 
compared  with  the  corresponding  values  of  pressure  calculated  fi'om  Eqs.  (2.3.18)  and  (2.3.28). 
To  this  end,  we  first  use  these  equations  to  calculate  values  of  the  pressure  along  the  isentrope 
passing  through  (v,,,  p=0) 

Because  shock  temperature  must  be  known  before  we  can  use  Eq.  (2.3.27),  we  first  use 
Eq.  (2.3.19)  to  evaluate  the  integral,  say  I,  on  the  right-handed-size  of  Eq.  (2.3.21),  and  thereby 
obtain  an  equation  for  calculating  along  the  Hugoniot  curve.  Because  r=2,  this  integral  I  can 
be  written  in  closed  form  as  the  sum  of  three  terms,  I  =  Ij+Ij+Ij,  with  I„  I2,  and  I3  expressed  in 
terms  of  V=v/vo  and  V=  (1-B)+BV,  as 

A  =wV2(1-D  (2.3.29.1) 

Equations  (2.3.21)  and  (2.3.29.1)  through  (2.3.29.3)  were  used  with  the  appropriate  values  of 

2(av,f  - 

I2  =  - ^[(K-l)-(2-5)lnF+(B-l)((T)  '-!)]  (2.3.29.2) 


4  =  (-^)'[0.5(1-F  )+(3-25)(F-l)+(3-5)(S-l)lnF-(5-l)2((F)-'-l)]  (2.3.29.3) 

Vo=0.3584  cmVg,  a=5.33  mm/psec,  and  B=1 .34  for  A1  to  calculate  values  of  Ty  in  the  reduced 
volume  range,  0.77  <  1 .  These  values  of  Th  and  the  corresponding  values  of  Pn  were  then 

used  in  Eq.  (2.3.27)  with  To=298K  to  calculate  values  of  the  pressure  Pj  along  the  (Vq,  p=0) 
isentrope.  The  energy  change  along  this  isentrope,  e^  -  eo,  was  also  calculated  using  Eq.  (2.3.6) 
with  r=2.  The  values  of  Tjj,  p^,  and  p^  and  Oj  -  eg,  calculated  in  this  way  are  listed  in  Table  2.1. 
The  following  equation  for  the  (vg,  p=0)  isentrope,  obtained  from  the  first  law  of 
thermodynamics  as 


Py  =  -^fv\^)dv 

J  dv 


(2.3.30) 


with 


giV)  =  , 

B  y  p 


(2.3.31) 

was  also  used  to  calculate  values  of  p^  at  the  same  values  of  V  to  provide  a  check  on  the  shock 
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temperature  calculations.  Both  of  these  sets  of  values  for  Pj  were  found  to  be  in  excellent 
agreement. 


The  equation  for  calculating  p^  from  the  Morse  potential  follows  from  Eqs.  (2.3.18)  and 
(2.3.28)  as 


Ps 


- ^[exp(26>')-exp(6y)]  ; 

3^2/3 


(2.3.32) 


the  equation  for  the  slope  of  this  isentrope  is  obtained  by  differentiating  Eq.  (2.3.32)  with  respect 
to  V,  yielding 


2  dp,  AbU.  2b^U. 

■■^0  "  - ^[exp(26y)-exp(6y)]  +  - — [2exp(26y)-exp(6y)]  (2.3.33) 


Recalling  that  the  isentrope  and  the  Hugoniot  curve  are  tangent  at  (vq,  p=0),  we  then 
equate  their  initial  slope  (dpH/dv)o=-(a/Vo)^  to  that  obtained  by  setting  y=0  in  Eq.  (2.3.33)  to 
obtain  the  following  equation  for  the  parameter  b: 

b"  =  9aV2Uo  (2.3.34) 

Because  the  listed  value  of  Uq  for  A1  at  298°K  is  78  kcal/mol,  and  a=5.33  mm/|j,sec,  we  are 
forced  to  use  the  value  of  b=3.25  when  calculating  Pj  from  Eqs.  (2.3.32).  The  values  of  p^  and  e^- 
eo  calculated  from  Eqs.  (2.3.32)  and  (2.3.18)  using  these  values  for  the  parameters  are  also  listed 
in  Table  2.1. 


We  now  compare  the  two  sets  of  values  for  p^  and  e^-eo  listed  in  Table  2. 1  to  see  if  the 
Morse  potential  is  thermodynamically  compatable  with  our  EOS  for  A1  based  on  shock  wave 
data.  The  excellent  agreement  between  these  two  sets  of  values  for  p^.  and  e^-eo  shows  that  the 
Morse  potential  can  be  used  for  A1  to  define  realistic  values  for  the  reference  energy  o^Uiy) 
along  the  reference  curve  Pr=  Ps(v)  =-d  U /dv  in  pressure  regions  of  interest.  In  this  case,  we  can 
write  Eq.  (2.3.2)  for  A1  with  r=2  and  the  function  g(v)  defined  by  the  Morse  potential  as 

e  =  pv/2  +  U  (v/vo)-PsV/2  (2.3.35) 

The  equation  for  p^  can  then  be  obtained  by  combining  Eqs.  (2.3.35)  and  (2.3.5)  as 

Ph  =  (p,V  +  2(eo-  t/(v/vo))po)/(2V-l),  (2.3.36) 

and  Eq.  (2.3.27)  with  p=p3  can  be  used  to  calculate  T^.  The  values  of  Pn  and  Th  calculated  from 
the  Morse  potential  in  this  way  are  also  listed  in  Table  2.1. 
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Table  2. 1 .  Values  of  the  shock  pressure  Ph,  the  shock  temperature  T^,  and  the  pressure  p^  and 
change  in  specific  internal  energy  e^-eo  along  the  (Vq,  p=0)  isentrope  calculated  for  different 
values  of  the  reduced  volume  V=v/Vo  from  an  EOS  for  A1  based  on  shock  wave  data  and  from  the 
Morse  potential  with  Uo=78  kcal/mol  and  b=3.25. 


EOS  Values 

Morse  Potential  Values 

V 

Th 

Cs-eo 

Ps 

Ps 

Ss-eo 

Th 

Ph 

K 

(nim/p)^ 

kbar 

kbar 

(mm/p)^ 

K 

kbar 

1.00 

0.0 

298.0 

0.0000 

0.0 

0.0 

0.0000 

298.0 

0.0 

0.95 

45.5 

332.2 

0.0391 

45.4 

45.3 

0.0399 

332.0 

45.4 

0.90 

105.7 

387.1 

0.1733 

104.7 

104.2 

0.1706 

387.1 

105.3 

0.85 

186.2 

495.5 

0.4243 

181.2 

180.9 

0.4234 

494.8 

185.9 

0.80 

295.8 

712.4 

0.8322 

279.9 

281.0 

0.8325 

718.3 

297.3 

0.77 

380.9 

933.6 

1.1713 

352.0 

355.6 

1.1725 

954.8 

385.9 

After  comparing  the  corresponding  sets  of  values  for  the  thermodynamic  variables  listed 
in  Table  2. 1 ,  we  conclude  that  the  Morse  potential  provides  a  realistic  description  of  A1  for  our 
detonation  calculations. 


2.3.2  An  Equation  of  State  For  AP  Based  on  Shock  Wave  Data. 

Because  AP  is  expected  to  decompose  completely  in  the  early  stages  of  the  detonation 
process,  we  use  Eq.  (2.3.6)  to  write  its  e=e(p,v)  EOS  as 

e=eH+(p-PH)v/r(v). 


In  this  case,  when  the  e=e(p,v)  EOS  is  rewritten  in  the  form  given  by  Eq.  (2.1.14)  as 


AP 


(Ah?) 


pv 


AP 


f  Up 


■  AP 


+  Vn 


AP 


Sap(-^ 


The  function  gAp(V)  is  defined  over  the  volume  spanned  by  the  Hugoniot  curve  passing  through 
the  initial  state  (v=Vo,  p=0,  e=eo).  A  method  of  extending  this  volume  range  is  presented  in 
Section  2.3.3. 


Our  construction  of  an  e=e(p,v)  EOS  for  AP  is  based  on  the  shock-wave  data  for  AP 
given  by  F.W.  Sandstrom,  P.A.  Persson,  and  B.  Olinger  (SPB)'',  and  its  thermodynamic  data 
given  in  the  LLNL  Explosives  Handbook'^.  To  be  more  specific,  SPB  give  the  following 
Hugoniots  for  solid  density  AP  (po=1.95  g/cm^),  U=2.84+1.85u  for  0.21  <  u  <  0.50  and 
U=2.90+l  .59u  for  0.66<u<2.10,  as  well  as  the  value  for  the  bulk  sound  velocity  Co=2.84 
mm/psec;  the  LLNL  Handbook  gives  the  estimated  values  for  the  specific  heat  at  constant 
pressure  Cp,  and  the  linear  expansion  coefficient  a=/'‘(d//3T)p  as  Cp=0.31  cal/g°C  and 
a«40xl0'V°C.  Substituting  these  values  for  Cp  and  Cq  into  the  thermodynamic  identity 
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r 


^0  dv 


c  V  dT^ 

p  p 


(2.3.37) 


and  the  value  v‘(0v/dT)p=3a=1.2xlO‘V'’C  gives  the  value  of  F  for  AP  as  F^p  =0.746.  Because 
shock  velocities  of  interest  are  in  the  6  niin/|isec  region,  and  because  of  uncertainties  in  the 
thermodynamic  data,  we  decided  to  base  our  EOS  for  AP  on  the  U=2.90+l  .59u  Hugoniot  curve 
and  the  value  Fap=0.75.  We  also  assumed  to  be  constant,  and  used  these  values  of  F  and 
v‘(dv/dT)p  with  T=298K  in  the  thermodynamic  identity 


(2.3.38) 


to  calculate  the  value  for  as  (Cv)ap=0.302  cal/g  k. 


Having  assigned  these  values  of  the  shock  and  thermodynamic  parameters  to  AP,  the 
combination  of  Eqs.  (2.3.5),  (2.1.14),  and  (2.3.19)  allows  us  to  write  the  g  function  for  AP  as 

g^,(V)  =  1.8333(p„)^/^)^[1-6,^-6,(4)^]  M  (2.3.39) 

B  V  V 

with  a=2.9  mm/psec,  B=1.59,  b,=0.8436,  b2=0.1564  and  F=((l-B)vo+Bv).  Equation  (2.3.2), 
with  Fap=0.75  and  gAp(V)  given  by  Eq.  (2.3.39)  can  now  be  used  to  calculate  the  isentrope 
through  the  initial  state  (Vo,p=0).  For  notational  convenience,  the  subscript  AP  will  be  mostly 
omitted  from  now  on. 


2.3.2.I.  The  Isentrope  Through  the  Initial  State 

The  differential  equation  for  the  isentropes  is  readily  derived  by  equating  the  expression 
de=-pdv,  obtained  from  the  first  laws  of  theormodynamics,  and  the  expression  for  de  obtained  by 
differentiating  Eq.  (2.3.2): 

vdp+(r+\)pdv  =  -V^dv  .  (2.3.40) 

dv 

Formally  integrating  Eq.  (2.3.40)  from  an  initial  condition  (v=Vj,  p=Pj)  gives  the  equation 

^  =p,.(v/^‘+F(v^g(v)-v,^g(v,,))+F'/.  (2.3.41) 

with  7.  =  Jv  ^~^g(v)dv ,  which  gives  the  equation  for  the  isentrope  passing  through  the  initial 

state  (Vi=Vo,  Pi=0)  as 

^  =  Fv^g(v)+F^7.  (2.3.41.1) 

with  -  Jv  ^■‘g(v)rfv ,  because  g(vo)=0.  For  convenience  in  evaluating  Iq,  we  set 

^0  ~  (2.3.42) 

with  a,=1.8333(a/B)2,  a2=1.5466(a/B)^  a3=0.2867(a/B)^  and 
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V  _j_  -JL  ^  -i. 

/j  =  J'v  '^dv  ,  I2  =  Jv  *V  ^dv  ,  -  jv  '* F  ^dv  . 


Our  construction  of  the  isentrope  through  the  initial  state  then  proceeds  as  follows. 
We  first  perform  the  itegration  in  I,  to  obtain  the  equation 

A  =  tK)''' 

3 


(2.3.43) 


with  V=Vo/v,  and  change  variables,  by  setting  v=X‘‘  and  F=B(XEBi)  with  Bi=(B-l)Vo/B.  With 
this  change  of  variables,  I2  and  I3  become 


4  /•  X^dX 


(2.3.44) 


and 


/  =  J_  f 


(2.3.45) 


We  then  use  the  identities 


-d- 


dX  ^  AX^dX 


XiX^^-B^)  X\X*-B^)  ix‘^-B^y 


and 


1  ,  1 


-] 


X\X^-B^)  iX^-B^) 

to  transform  the  integral  I3  in  two  steps,  the  first  by  parts  and  the  second  by  using  partial 
fractions,  into  the  equation 


(2.3.46) 


(2.3.47) 


3/4 


1]  - 


(5-l)Vo''  (5F-(5-1))  4(5 -1)v 


(2.3.45)' 


At  this  stage,  it  is  convenient  to  rewrite  the  equation  for  Iq  as 


3  *  “  5-1  5F-(5-l) 


-1]  +  ( 


4(5-1) 


■«2)''o4  • 


(2.3.48) 


We  finally  use  partial  fractions  to  perform  the  integration  in  Ij  and  obtain  the  following  equation: 
h  =5-‘(._^..  1  In-^ - ?—  +  In - +  2(tan-‘(^)^  -  (2.3.44)' 


(5 -IK 


1  1 

B  B 


which  together  with  Eq.  (2.3.45)'  defines  the  integral  I3  as  a  function  of  the  reduced  volume 
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V=v/Vo. 


Equation  (2.3.39)  and  Eqs.  (2.3.44)'  and  (2.3.48)  were  used  to  calculate  values  of  gAp(V), 
I2  and  Io(V)  over  the  reduced  volume  range  0.62^V<  1,  and  then  Eq.  (2.3.41 .1)  was  used  with 
these  values  of  gAp(V),  Io(V),  and  F ^=0.75  to  calculate  values  of  the  pressure  on  the  isentrope 
over  the  same  volume  range.  These  values  of  gAp(V),  Io(V),  and  Pj  are  given  in  Table  2.2,  as  well 
as  the  values  of  the  temperature  along  the  isentrope  calculated  from  the  equation 

T  =  r.(^)®”  (2.3.49) 

V 

with  Vi=Vo  and  Ti=To. 


To  complete  this  series  of  EOS  calculations  for  AP,  we  calculated  values  of  the  pressure 
Ph  and  temperature  Tjj  along  the  Hugoniot  in  the  (V-p)  plane  over  the  same  reduced  volume 
range,  0.62<  V <  1 .  These  values  of  Pn  and  Th  are  also  given  in  Table  2.2.  The  values  of  pn  were 
calculated  from  Eq.  (2.3.19)  with  Vo=0.513  cmVg,  a=2.90  mm/psec  and  B=1.95,  and  the  values 
of  Th  were  calculated  by  substituting  the  corresponding  values  of  T^,  p^,  and  Pn  into  the  equation 


■  H 


T  +P'£^-Zi 

^  Pore, 


(2.3.50) 


with  po=1.95  g/cm^,  rAp=0.75  and  (CJap^I  •264x10'^  kbar  cmVg  K. 

Examination  of  the  values  of  Th  in  Table  2.2  leads  us  to  question  this  treatment  of  AP 
because  we  do  not  expect  the  shock  temperature  to  reach  such  high  values  in  the  300-400  kbar 
region.  The  reason  for  these  high  temperatures  is  clearly  seen  hy  considering  Eq.  (2.3.50).  As  V 
decreases  and  (pn-Ps)  increases,  the  difference  (Th-T^)  becomes  magnified  because  Cy  and  F  are 
assumed  to  be  constant  and  F<1.  We  suggest  that  a  task  to  investigate  the  validity  of  these 
assumptions  be  made  a  task  in  any  additional  work  that  may  be  performed  on  explosive 
compositions  containing  AP. 


2.3.2.2.  The  Family  of  Isentropes  Intersecting  the  Hugoniot  Curve. 

Because  of  our  assumption  about  the  behavior  of  reacting  AP  behind  a  shock,  we  must 
have  a  method  for  calculating  the  relationship  between  the  pressure  and  volume  along  isentropes 
that  intersect  the  shock.  An  equation  providing  this  relationship  along  the  isentrope  passing 
through  a  point  (Vh,  Ph,  Th)  on  the  Hugoniot  is  readily  obtained  by  arranging  Eq.  (3.2.50),  with 
Ts=Th(Vh/Vo)°’^  as 

FC  T  V 

Ps  =  Ph-9,-^  (1  (2.3.50)' 


with  Vh=Vh/vo. 


Table  2.2.  Values  of  v^gAp(V),  I(/(vo)ap  the  pressure  p^  and  temperature  Tj  along  the  (v^,  p=0) 
isentrope;  the  shock  pressure  Pn  and  the  shock  temperature  Th  for  different  values  of  the  reduced 
volume  V=v/Vo. 
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V 

-vfgAp(v) 

(mm7jj,sec)2 

V('^o)aV 

(mm/psec)^ 

Ps 

kbar 

T 

K 

Ph 

kbar 

Th 

K 

1.0 

0.000 

0.0 

298.0 

0.0 

298.0 

0.95 

0.616 

9.68 

309.7 

9.68 

309.7 

0.90 

1.368 

23.09 

322.5 

23.19 

327.1 

0.85 

2.302 

41.95 

336.6 

42.42 

358.2 

0.80 

3.496 

68.95 

352.3 

70.52 

420.2 

0.75 

5.068 

108.57 

369.8 

112.94 

547.1 

0.70 

7.226 

0.866 

168.68 

389.4 

179.86 

812.9 

0.65 

10.351 

1.347 

264.30 

411.6 

291.82 

1379.6 

0.62 

12.990 

1.736 

350.34 

426.5 

397.80 

2018.8 

2.3.3.  A  Method  to  Extend  Equations  of  State  Based  On  Shock  Wave  Data 

We  describe  here  another  set  of  E,P,V  and  T,P,V  equations  of  state  of  AP,  which  differ 
from  those  derived  in  Section  2.3.2  in  the  following  ways.  With  respect  to  the  E,P,V  equation  of 
state,  the  two  are  identical  for  the  region  V  <  Vq;  in  this  section,  we  extend  the  equation  of  state  to 
the  region  V>Vo,  which  was  not  done  in  Section  2.3.2.  The  T,P,V  equations  of  state  are  different 
in  that  Cy=constant  is  assumed  in  Section  2.3.2,  whereas  in  this  section,  we  use,  instead,  the 
experimental  coefficient  of  volume  expansion  and  specific  heat  at  constant  pressure  along  the  1 
bar  isobar.  The  derivations  provided  below  closely  follow  the  earlier  theoretical  results  of  Enig.^ 

Enig®  used  the  following  Mie-Gruneisen  form  of  the  E,P,V  equation  of  state,  which  we 
now  apply  to  ammonium  perchlorate: 

P-P,(V)=(yA^)[E-E,(V)],  (2.3.51) 

where  P,V,  and  E  are,  respectively,  the  pressure,  specific  volume,  and  specific  internal  energy,  y 
is  the  Gruneisen  coefficient  (assumed  to  be  a  constant  for  all  P  and  V),  and  subscript  r  refers  to 
the  reference  state  to  be  specified  below.  It  should  be  noted  that  Eq.(2.3.5 1)  is  a  first-order 
Taylor  expansion  in  the  neighborhood  of  the  reference  curve,'  and 

Y=V(0P/dE)v=constant.  (2.3.52) 


In  order  to  evaluate  y,  we  note  that  the  definition  of  the  sound  speed  c, 

c2=v'[P(ap/aE)v-(ap/av)E],  (2.3.53) 

may  be  rewritten  as 

cW=[p+(aE/0V)p](ap/aE)v 
=(Cp/pv)(ap/aE)v, 

where 

p=(W)(aV/aT)p 


(2.3.54) 

(2.3.55) 
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and  Cp  are,  respectively,  the  coefficient  of  thermal  expansion  and  specific  heat  at  constant 
pressure.  From  Eqs.  (2.3.52)  and  (2.3.54),  it  follows  that 

Y=Pc'/Cp,  (2.3.56) 

where  P,  c,  and  Cp  may  be  chosen  at  the  initial  state  (Po,Vo,Eo,Tp),  i.e., 

Y=Yo=PoCo2/Cp,o.  (2.3.57) 

and  the  subscript  o  refers  to  the  initial  state  Po=l  bar,  To=298°K. 

Thus  far,  the  reference  state  has  not  been  specified.  Usually,  for  V<  V^,  the  Hugoniot 
originating  at  the  initial  state  Ep,Pp,Vo  is  chosen  as  the  reference  state;  we  also  shall  do  that.  For 
the  usual  application  of  an  E,P,V  equation  at  high  pressures  as  encountered,  e.g.,  in  the  shocking 
of  a  metal,  the  thermodynamic  states  are  limited  to  W<Y„.  However,  when  the  shocked  state  is 
subsequently  expanded  isentropically  by  a  rarefaction  wave  originating  at  a  free  surface,  the  end 
state  may  have  V>Vo.  If  P^V)  and  E^V)  refer  to  the  Hugoniot,  V  is  not  defined  for  V>Vo.  The 
usual  way  this  dilemma  is  treated  is  formally  to  extend  the  P,V  Hugoniot  relationship  to  V>V„, 
i.e.,  into  the  "rarefaction  shock"  region.  It  is  argued  that  the  second-order  tangency  of  the 
Hugoniot  and  the  isentrope  at  (Po,Vo,Eo)  justifies  this  extrapolation  if  the  volume  increase  in  the 
expanded  state,  produced  by  the  irreversible  shock  heating,  is  small.  This  argument  may  be 
acceptable  for  the  calculation  of  material  motion,  but  probably  would  lead  to  large  errors  in  the 
calculation  of  the  temperature  of  the  shocked  state.  It  will  be  shown  below  that  the  temperature 
of  the  shocked  state  is  a  sensitive  fimction  of  the  T,V  variation  along  P=Po  and  the  slope  of  the 
P,V  isentropes  at  P=Po.  For  this  reason,  we  shall  use  the  isobar  P=Po  as  the  pressure  reference 
state  Pj  and  Ej(V)=E(Po,V)  as  the  internal  energy  reference  state. 

2.3.3.I.  Region  1:  V 

In  the  low-pressure  region  V^Vp,  Enig  defined^ 

Pr(V)=Po  (2.3.58) 

E,(V)=E(P„,V).  (2.3.59) 


By  integration  of  the  thermodynamic  identity 

(aE/aV)p=Cp/pV-P  (2.3.60) 

along  P=Po,  we  find 

E(P„,V)=E„+/vjCp/pV’-PJdV’,  (2.3.61) 

where  Eo=E(Po,Vo),  the  specific  internal  energy  at  the  initial  state,  is  an  arbitrary  constant  that  can 
be  set  equal  to  zero  if  desired.  If  Cp=Cpo  and  P=P„  on  P=Po,  Eqs.  (2.3.55)  and  (2.3.61)  yield 

T(P„,V)=T„+(1/PXVA^„)  (2.3.62) 
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EXV)=E(P„,V)=E„+(Cpypjf«(VA^J-P„(V-V„^  (2.3.63) 

It  should  be  noted  that  though  all  terms  involving  P^  are  negligible,  they  will  be  kept  in  order  that 
all  thermodynamic  relationships  to  be  derived  are  formally  exact.  Equation  (2.3.51)  becomes 

E=E„-P,(V-VJ+(Cp,yp„)««(VA^„)+(P-PJV/Y  for  V„.  (2.3.64) 

Enig*'*  has  shown  that  the  thermodynamic  identity 

[p+(aE/av)p](aT/0P)v-(aE/ap)v(aT/av)p=T  (2.3 .65) 

yields  the  characteristic  equations 

dP/dv-i-[p-i-(aE/av)p]/(aE/ap)v=o  (2.3 .66) 

dT/dV+T/(aE/aP)v=0,  (2.3.67) 


which  with  an  E,P,V  equation  of  state  and  a  nonisentropic  initial  curve  along  which  the  tempe¬ 
rature  is  prescribed,  e.g.,  Eq.  (2.3.62),  leads  to  a  T,P,V  equation  of  state  as  described  below. 
The  solution  of  Eq.  (2.3.66)  yields  the  isentropes 


C,=[(P-Po)V+CpypJV\  (2.3.68) 

where  the  constant  of  integration  C,  is  different  for  each  isentrope,  i.e.,  C,  varies  with  the  entropy 
S.  The  solution  of  Eq.  (2.3.67)  along  an  isentrope  is 

C2=TY\  (2.3.69) 

where  Cj,  the  constant  of  integration,  varies  with  S.  Substitutions  of  P=Pp  into  Eq.  (2.3.68)  and 
Eq.  (2.3.62)  into  Eq.  (2.3.69)  yield,  respectively,  C,  and  C2  as  functions  of  the  parameter  V  (the 
specific  volume  along  P=Po): 

Ci=Cp,oVVp„  onP=P„  (2.3.70) 

Q=[T„+(l/p.)(»(VWJ]V’'  onP=P<,.  (2.3.71) 

For  convenience,  replace  V  in  the  above  equations  by  z.  Equations  (2.3.68)  and  (2.3.69)  are  now 

[(P-PJV+Cp,ypJV^=Cp.„zVp„  (2.3.72) 

TVv=[T„+(l/pXzA^„)]zV  (2.3.73) 

Elimination  of  the  parameter  z  between  Eqs.  (2.3.72)  and  (2.3.73)  yields  the  P,V,T  equation  of 
state: 


T(P,V)=x[T.+(l/p.)te(x''’VAt.)],  (2.3.74) 

where 
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x(P,V)=l+P„(P-P„)V/Cp,„. 


(2.3.75) 


In  particular,  on  V=Vo: 

T(P,VJ=TJl+p„(P-PJVyCpJ{l+(l/YpoT„)^«[l+Po(P-Po)Vo/Cp.o]}  (2.3.76) 

23.3.2.  Region  2:  V<V„ 

For  V<Vo,  Enig’  chose  the  Hugoniot  pressure  Ph(V)  as  the  reference  pressure  P^CV), 
where  the  origin  of  the  Hugoniot  is  the  initial  state  (P^jV^jEJ.  Thus,  Eq.  (2.3.51)  is  now 


P=Ph(V)+(yA^)[E-Eh(V)].  (2.3.77) 

Substitution  of  Eq.  (2.3.77)  into  the  isentropic  relationship 
dE+PdV=0  (2.3.78) 

and  use  of  the  shock  condition 


Eh-E=(Ph+Po)(Vo-V)/2,  (2.3.79) 

yields 

dE/dV+YEA^=Y[E„+‘/2Po(V„-V)]m‘/2YPH(Vo-V)A^-PH,  (2.3.80) 

whose  general  integral  is 


C3=(VA^o)lE-E„-»/2P„V„+y2YPoV/(Y+l)]+E,+y2P<.Vy(Y+l)+I(V),  (2.3.81) 

where 

I(V)=i/2V„-vj^<,vV(2+Y-YVA)PH(v)dv  (2.3.82) 


and  C3  is  the  constant  of  integration  defining  the  isentrope  in  the  E,V  plane.  In  terms  of  the  as 
yet  unknown  Hugoniot  Ph(V),  the  isentropes  in  the  P,V  plane  are  obtained  by  substitution  of  E 
from  Eq.  (2.3.77)  into  Eq.  (2.3.81):  , 


V 


P  V 


P.F_ 


C'3=(— )''[(P-P^)-+V4(P^+P„)(F -E)-'/2P„F  +*/2Y-^]+£„+'/2^+/(F). 


Y+1 


Y+1 


(2.3.83) 


To  this  point,  the  analysis  for  Region  2  is  general  in  terms  of  an  as  yet  nonprescribed 
Hugoniot  curve  PyfV).  For  simplicity  of  analysis,  Enig^  limited  himself  to  those  unreacted 
explosives  (as  well  as  metals)  that  satisfy  the  experimental  shock  condition 


U=Co+au,  (2.3.84) 

where  U  and  u  are,  respectively,  the  shock  velocity  and  particle  velocity,  and  c,,  and  a  are 
empirical  constants. In  theory,  the  c^  of  Eq.  (2.3.84)  should  be  the  same  as  the  initial  sound 
speed  of  Eq.  (2.3.57),  but  in  practice  it  is  not  because  Eq.  (2.3.84)  generally  does  not  hold  as  u-0. 
It  follows  from  the  Rankine-Hugoniot  conditions  that 


PH(V)=P„+c„^(V„-V)/[V„-a(V„-V)]l 


(2.3.85) 
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Substitution  of  Eq.  (2.3.85)  into  Eq.  (2.3.82)  gives 

I(V)=/2P„V„[(Y+2)Rv^V(Y+l)-R^KY+l)‘']+c^A^''(Y+l)(a-l)-^J^^^  (2.3.86) 

where 

R=VA^„  (2.3.87) 

J(V)=jr^-'^(mt-l+n/t)(l-t)-Mt  (2.3.88) 

m=!/2Ya/[(a-l)(Y+l)]  (2.3.89) 

n=‘/2(Y+2)(a-l)/(Y+l)a  (2.3.90) 

A=l-a-‘<l  (2.3.91) 

B(V)=AVyV<l,  (2.3.92) 


and  the  transformation  t=(l-a'*)Vyv  has  been  used.  It  should  be  noted  B(V)^  1,  where  the 
equality  holds  only  as  V/Vp-*  l-a'*  as  Ph-^°°  in  Eq.  (2.3.85). 

Using  partial  fractions,  Eq.  (2.3.88)  becomes 

J(V)=JiWt-^[A,t-'+A2(l-t)-'+A3(l-t)-2]dt,  (2.3.93) 

where 


Ai=n,  A2=n-m,  A3=m+n-l. 


(2.3.94) 


But 


jB(v)t-Y(i.t)-2dt=B'-Y(l-B)-'-A‘-Y(l-A)-'+Yjr^f^(l-t)-'dt,  (2.3.95) 


which  when  substituted  into  Eq.  (2.3.93)  yields 

J(V)=[A,/Y-AA3/(l-A)]A-^-[A,/Y-BA3/(l-B)]B-^+(A2+YA3)|®(^>r(l-t)-'dt.  (2.3.96) 

By  Eqs.  (2.3.83),  (2.3.86),  and  (2.3.96),  the  evaluation  of  the  constant  C3,  which 
characterizes  an  isentrope,  entails  the  evaluation  of  the  integral  in  Eq.  (2.3.96).  To  evaluate  the 
integral,  we  note  that 


J®(^t-^(l-t)-Mt=J«<^^-Mt/[t^(l-t)]=Ji(''H^-MtEt-'+jrH^-^(l-t)-'d^ 


= jB(v)tN-Y(  j  _^).idt+53(B^-'^-T'^'-AN-''-T^+>)/(N-k-Y+l ) 

=J®(^¥-'(l-t)-’dt+E(B^-''-A^-'')/(s-k),  (2.3.97) 

where  the  summation  is  taken  from  k=l  to  k=N,  N  is  the  smallest  non-negative  integer  greater 
than  Y-1, 


30 


ENIG  TR  95-1 


s=N-y+1>0, 
and  the  integral*^ 


(2.3.98) 


jrH"'(l-t)-'dt=s-'B^2F,(l,s;l+s;B)-s-'A^2F,(l,s;l+s;A),  (2.3.99) 

is  given  in  terms  of  Gauss's  hypergeometrie  series.’’  Note  that  if  N=0,  the  summation  term  on 
the  right-hand  side  of  Eq.  (2.3.97)  is  set  equal  to  zero.  The  exaet  derivation  of  the  P,V  isentrope 
characterized  by  some  value  of  Cj  (a  fimction  of  the  entropy)  is  complete. 

In  order  to  avoid  the  need  to  evaluate  Gauss's  hypergeometrie  series,  the  integral  in  Eq. 
(2.3.96)  can  be  accurately  approximated  by  fitting  (1-t)''  by  a  polynomial  in  t  over  the  range 
0<l-a''=A<t<B(V)=AVoA^<l  and  carrying  out  the  resulting  simple  integration.  A  single 
polynomial  fit  can  be  made  to  work  for  almost  all  explosives  over  a  single  interval  A^i„<t<Bn,a,j, 
where  A„,i„-0.2  and  B„^-0.9: 

(1-t)-'  = 

where  ao  =  13.29242189,  a,  =  247.2782554,  a2  =  1745.42271 1,  aj  =  6618.403805,  = 

14505.10394,  aj  =  18455.64985,  a^  -  -12662.06958,  and  a,  =  3633.143062. 

The  solution  of  Eq.  (2.3.67)  along  an  isentrope  is 

C4=TV^  (2.3.100) 

where  the  subscript  4  is  used  in  this  region  to  distinguish  the  isentrope  from  that  in  Region  1 
given  by  Eq.  (2.3.69).  Substitution  of  Eq.  (2.3.76)  into  Eq.  (2.3.100)  yields  C4  as  a  function  of 
the  parameter  P  on  the  non-characteristic  (non-isentropic)  line  V=V„;  then,  for  convenience,  P  is 
replaced  by  the  new  parameter  y.  By  Eq.  (2.3.100)  the  temperature  along  the  isentrope  is  now 
given  by 

TV^=T„V/[l+p„(y-PJVycJ{l+(l/YPoT„)««[l+Po(y-Po)Vo/^^^^^  (2.3.101) 

Substituting  V=Vg  into  Eq.  (2.3.83)  yields 


C3=E„+(P-P„)VA  onV=V„.  (2.3.102) 

Replacing  P  by  y  in  Eq.  (2.3.102)  and  substituting  the  result  into  Eq.  (2.3.83)  gives  the 
P,V  isentrope: 

o  ^  o 

Equations  (2.3.101)  and  (2.3.103)  define  the  T,P,V  equation  of  state  for  Region  2  in  terms 
of  the  parameter  y.  Thus,  for  known  P  and  V,  y  is  found  from  Eq.  (2.3.103)  (remembering  that 
Ph  is  given  by  Eq.  (2.3.85))  and  then  substituted  into  Eq.  (2.3.101)  to  yield  T. 

2.3.3.3.  Connecting  the  Isentropes  of  Regions  1  and  2  at  V=V„ : 


— +'/2(P„+P  )(V  -V)-V2P  V  +V2yp  — — 

Y  ^  H  o  J  00  I 


P  V 


(2.3.103) 
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The  P,V  isentropes  in  Regions  1  and  2,  which  are  given,  respectively,  by  Eqs.  (2.3.68) 
and  (2.3.83),  meet  at  V=V„.  Substituting  in  Eqs.  (2.3.68)  and  (2.3.73)  and  eliminating  P 
between  the  resulting  equations  give 


C,=[Y(C3-E„)+CpypJVJ. 


(2.3.103.1) 


2.3.4.  A  Thermodynamic  Model  For  AI2O3. 


An  examination  of  the  shock  wave  data  for  corundum’^,  a-AljOj,  with  a  density  of 
(Po)i=3.98  g/cm^,  and  our  shock  pressure  range  of  interest  led  us  to  conclude  that  an  adequate 
thermodynamic  treatment  of  a-AljOj  and  Y-AI2O3  can  be  based  on  the  assumption  that  the 
Hugoniot  curves  are  straight  lines.  In  this  case,  the  equations  for  the  pressure  and  energy  along 
these  Hugoniot  curves  can  be  written  as 


and 


P  =  (PoS).  K"'^)/  ,  *=1  O'-  2  , 


^  (Po^oX  f  .^2  . 

=  — ;; — ,  '=1  or  2  , 


(2.3.104) 

(2.3.105) 


where  (Cq);  denotes  the  constant  shock  velocities,  (e^).  =  (Aft^“)^,  ^  ,  and  the  Hugoniot  curves  are 
also  isentropes  because  Eq.  (2.3.105)  satisfy  the  equation  dej—pav-.  The  combination  of  Eqs. 
(2.3.104),  (2.3.105)  and  the  e=e(p,v)  EOS  given  by  Eq.  (2.1.16)  gives  the  following  equations  for 
the  gi(V)  fimctions  for  AI2O3: 


g,(V)  =  pi, 


(S)? 

2 


2(rt-i)r ,  (r+2)F^ 
t  r  r 


/or  F^l  ,  /  =  1  or  2  . 


(2.3.106) 


To  proceed  with  our  thermodynamic  model,  we  assumed  that  a-Al203  and  Y-AI2O3  have  constant 
values  of  r  and  Cy.  With  literature  values  of  (Cp),=0.775xl0'^  kbar  cmVg  K, 

((5v/3T)p),=5. 108x10'^  kbar  cm^/g  K  for  a-Al203,  and  a  value  (Co),=9.64  mm/psec  suggested  by 
the  shock-wave  data,  Eqs.  (2.3.37)  and  (2.3.38)  give  the  values  r',=2.44  and  (Q), =0.764x10'^ 
kbar  cmVg  K,  for  a-Al203. 


Now  with  the  limited  amount  of  thermodyanamic  data  available,  we  have  to  formulate  a 
procedure  for  deciding  whether  a-Al203  or  Y-AI2O3  is  present  in  our  reacting  explosive 
compositions  containing  Al.  We  base  this  formulation  on  the  fact  that  a-Al203  is  the  most  stable 
form  of  AI2O3  in  the  standard  states  (p=0,  T=298,  v=(vo)i),  and  on  the  premise  that  the  Hugoniot 
curves  of  a-Al203  and  Y-AI2O3  intersect  at  a  point  (vj,  pj)  in  the  (v,p)  plane.  Then,  AI2O3  will 
be  assumed  to  be  in  the  a-form  when  p<Pj,  but  in  the  y-form  when  p>pj.  Our  justification,  for 
this  assumption  as  well  as  our  procedure  for  evaluating  unknown  thermodynamic  parameters  for 
Y-AI2O3  is  based  on  the  following  calculations  of  the  Gibbs  free  energies  Gj  and  G2  for  a-Al203 
and  Y-AI2O3 ,  respectively. 


Our  free  energy  calculations  are  based  on  the  thermodynamic  data  for  a-Al203  and  y- 
AI2O3,  which  are  listed  in  the  66th  Edition  of  the  Handbook  of  Chemistry  and  Physics  (HCP)**. 
The  values  given  for  the  heat  of  formation  and  entropy  for  a-Al203  in  the  standard  state  (p=l 
atm,  T=298.15  K)  are  (AH®), =-400.5  kcal/mol  and  S°=12.17  cal/deg  mol.  In  constrast  to  a-Al203 
,  a  range  of  values  of  (po)2  and  (AHf)2  are  given  for  Y-AI2O3,  and  no  value  for  S®  is  given.  To  be 
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more  specific,  as  Y-AI2O3  changes  from  its  amorphous  to  crystalline  form  in  the  standard  state, 
(po)2  varies  from  3.5-3.9  g/cm^  and  (AHf)2  varies  from  -390  to  395  kcal/mol.  Additional 
background  information  for  our  thermodynamic  model  of  AI2O3  was  the  following  experimental 
observation  from  a  paper  by  A.  J.  Brock  and  M.  J.  Pryor,'^:  the  only  oxidation  product  formed  by 
pure  A1  subjected  to  low  pressures  of  O2  in  the  300-425‘’C  temperature  range  is  amorphous  y- 
AI2O3.  This  observation  suggests  that  the  Y-AI2O3  formed  during  the  detonation  process  will  be 
in  a  somewhat  amorphous  state;  it  allows  us  to  make  an  estimate  of  the  standard  entropy  S®  of 
this  amorphous  Y-AI2O3. 

Our  estimate  of  S®  is  based  on  the  assumptions  that  a-Al203  and  amorphous  Y-AI2O3  have 
the  same  values  of  Cp,  (Cp)i=(Cp)2,  and  that  Gi=G2  at  300°C.  In  this  case,  the  equation 

-  TSI  (2.3.107) 

with  T=573  K,  the  values  for  (AHf),  and  8°  given  above,  and  a  value  of  (AHf)2~392  kcal/mol 
gives  a  value  for  the  standard  entropy  of  amorphous  Y-AI2O3  of  82=27.00  cal/deg  mol. 

At  this  stage  in  our  treatment  of  AI2O3,  we  have  two  adjustable  parameters,  pj  and  (pojj, 
whose  values  must  be  specified  before  this  model  can  be  utilized.  Because  values  for  these 
parameters  are  not  known  at  this  time,  we  first  present  a  procedure  for  completing  this  model 
when  Pj  and  (po)2  are  prescribed,  and  then  exemplify  it  by  using  particular  values  for  pj  and  (po)2. 

Our  procedure  for  closing  our  thermodynamic  treatment  of  AI2O3,  when  pj  and  (po)2  are 
prescribed,  is  as  follows.  First,  use  the  value  of  pj  in  Eq.  (2.3.104)  with  i=l  to  determine  the 
value  of  Vj,  and  use  this  value  of  Vj  in  the  following  equation  with  i=l, 

(9.  =  T„(^)'’'  ,  (2.3.108) 

to  calculate  the  value  of  (Tj)„  at  (vj,  pj).  Next,  use  the  condition  G,=G2  at  the  point  (vj,  pj)  to 
calculate  (Tj)2,  and  use  Eq.  (2.3.108)  with  j=2  to  calculate  the  corresponding  value  of  r2. 

Because  (Co)2  is  known  from  Eq.  (2.3.104)  when  pj  and  (Vo)2  are  specified,  and  (Cp),=(Cp)2  by 
hypothesis,  Eq.  (2.3.37)  can  be  used  to  calculate  the  corresponding  value  of  ((6v/5T)p)2  for  y- 
AI2O3.  We  will  now  carry  out  this  procedure  in  detail  for  the  case  when  pj=100  kbar  and 
(Po)2=3.52  g/cm^ 

Thus,  we  first  set  Pj=100  kbar  in  Eq.  (2.3.104)  with  (poCg), =3.698x10^  kbar  to  obtain  the 
corresponding  value  of  Vj=0.2445  cmVg.  Equation  (2.3.108)  with  i=l  and  ri=2.44  then  gives  a 
value  of  (Tj)i=318.6  K  and  Eq.  (2.3.104)  with  i=2  and  (Vo)2=0.2841  cmVg  gives  the  value  of 
(Co)2=4.515  mm/|isec.  We  next  write  the  equations  for  G,  and  G2  at  (vj,  pj)  as 

G,  =  (Aff/),  +  (2.3.109.1) 

and  use  the  equation  obtained  by  setting  Gi=G2  to  calculate  a  value  of  (Tj)2=613  K,  which  is  used 
in  Eq.  (2.3.108)  with  i=2  to  calculate  a  value  of  r2=4.80.  Because  the  value  of  (Tj)2=340°C  lies 

G,  =  (Ah/),  *  (2.3.109.2) 
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towards  the  bottom  of  the  temperature  range  reported  for  the  formation  of  amorphous  y-AljOj  at 
low  pressures,  it  is  reasonable  to  choose  Pj=100  kbar  and  (po)2=3.52  cmVg  as  the  values  for  these 
parameters  in  a  first-generation  thermodynamic  model  for  AI2O3  developed  in  T^lsk  1.  The 
calculated  value  of  r2=4.80  leads  to  the  conclusion  that  in  this  treatment  of  AI2O3  the  thermal 
coefficient  of  amorphous  Y-AI2O3,  ((dv/5T)p)2,  is  ten  times  greater  than  the  thermal  coefficient  of 
a-Al203,  ((dy/dT\\. 


In  any  future  work  to  be  performed  on  explosive  compositions  contains  Al,  the  method 
presented  in  the  previous  section  should  be  used  to  extend  these  EOS  for  a-Al203  and  amorphous 
Y-AI2O3  beyond  the  volume  ranges  spanned  by  their  Hugoniot  curves.  Melting  of  AI2O3  should 
also  be  considered. 

2.3.5.  Polytropic  Equations  of  State  for  the  Detonation  Products 

The  limitations  of  the  polytropic  EOS  led  to  the  development  and  implementation  of 
more  realistic  EOS  for  detonation  products^®’^'.  More  basic  EOS  for  detonation  products,  based 
on  the  work  of  J.D.  Weeks,  D.  Chandler,  and  H.C.  Anderson^^  have  been  developed  by  F.H. 

Ree^^  and  W.B.  Brown^'’ ,  and  for  detonation  products  containing  Al  by  H.D.  Jones  and  F.J. 
ZerillP  .  These  EOS  address  both  the  spherically  averaged  and  multipole  components  of  the 
temperature-dependent  pair  potential  between  the  molecular  species.  Use  of  these  EOS  was 
considered,  but  we  decided  that  there  was  insufficient  time  to  adopt  them  in  Phase  I.  For  this 
reason,  the  polytropic  EOS  was  considered  to  be  adequate  for  our  detonation  modelling  in  Phase 
I.  Consequently,  the  polytropic  E=E(p,Vf,mH2o)  EOS  of  our  AP/Al  mixture  has  already  been 
introduced  in  Section  2.1,  and  the  polytropic  E=E(p,Vf,mij2o,  inco2,  inc2H6)  EOS  for  our 
CaH^NpOyAl  mixture  can  be  seen  in  Eq.  (2.2.20).  These  polytropic  EOS  extend  those  usually 
encountered  in  the  literature  by  accounting  quantitatively  for  the  reactions  of  the  detonation 
products  with  Al.  The  reaction  coordinates  used  to  describe  the  changing  masses  of  some  of  the 
fluid  constituents  in  the  detonation  products,  however,  were  only  included  in  the  standard  states 
term  containing  heats  of  formation  as  shown  in  Eq.  (2.2.20). 

Because  these  E=E(p,Vf,mH2o)  and  E=E(p,Vf,mH2o,  nico25  JW:2h6)  EOS  are  incomplete  and 
provide  no  information  about  temperature,  they  are  only  adequate  for  treating  the  detonation 
process  when  the  temperature  need  not  be  accounted  for  explicitly.  Consequently,  they  are 
adequate  for  a  detonation  model  when  the  solid  components  behind  the  shock  are  assumed  to  be 
constrained  to  their  own  isentropes,  but  not  when  a  solid  component  such  as  AI2O3  is  assumed  to 
attain  thermal  equilibrium  with  the  detonation  products.  Thus,  it  was  necessary  to  construct  a 
T=T(p,Vj)  EOS  for  detonation  products  to  accomodate  our  assumption  that  Tf=Ti.  This 
construction  proceeds  as  follows. 

2.3.5.I.  A  p=p(T,V)  EOS  of  Detonation  Products 

To  be  consistent  with  the  polytropic  EOS  introduced  in  previous  sections,  we  write  the 
complete  E=E(Sf,Vf,  ^j)  EOS  for  detonation  products  as 

E  =  ,  (2.3.110) 

so  that  the  extents  of  reaction  appear  only  in  the  standard  state  term  Ec(^i).  Now,  because  we 
are  dealing  with  polytropic  media  only,  it  is  convenient  to  omit  the  subscript  f  In  addition. 
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where  it  is  convenient  for  our  purpose,  we  will  only  consider  one  extent  of  reaction  coordinate 


In  this  case,  the  thermodynamic  identities  for  T  and  p  can  be  written  as  T=0Ei/5S  and 
p=-9Ei/3V,  the  EOS  for  T  and  p  in  terms  of  S  can  be  written  as  T=T(S,V)  and  p=p(S,V),  and  S 
can  be  eliminated  between  these  EOS  to  show  that  T=T(p,V).  Thus,  when  E=E(S,V,^)  is  written 
as  E=Ec(0+Ei(S,V),  the  partial  derivatives  (0T/0Op,v’  (^T/dO  s.v>  s,v,  (5S/dO  x.v  are 

all  zero,  and  it  follows  from  the  following  identities  for  the  chemical  affinity  A“, 


-A 


dE  _  dE  dE 


(^) 

dp 


(2.3.111) 


that  -A=(dE/dO  p  v, .  Consequently,  the  ^  terms  in  Ec(Q  are  only  associated  with  the  heats  of 
formation  of  the  species  taking  part  in  the  reaction. 


To  proceed  with  our  construction,  we  consider  the  case  when  E,  is  a  separable  function  of 
S  and  V  given  by  the  equation 

E^  =  V-"E^(S)  ,  (2.3.112) 

where  n  is  a  constant  parameter,  and  we  differentiate  this  expression  for  Ej  to  obtain  equations 
for  p  and  T.  The  identity  for  p  gives  the  equation 

—  =  >  (2.3.113) 

n 

which  together  with  Eq.  (2.3.1 10)  allows  us  to  write  the  E=E(p,V,^j)  EOS  as 

E  =  E^a)  +  ^  (2.3.114) 

n 


and  identify  the  parameter  n  as  n=Y-l .  The  identity  for  T  gives  the  equation 

T  =  V'"dE^/dS  , 

which  leads  to  the  following  equation  for  C^, 

dEJdS 

Cy  =  - - -  . 

d^E^/dS^ 


(2.3.115) 

(2.3.116) 


We  next  combine  Eq.  (2.3.1 15)  with  the  equation  obtained  by  differentiating  Eq.  (2.3.1 13) 
partially  with  respect  to  T,  to  obtain  the  following  equation: 

.  (2.3.117) 

Because  (dS/dOjv“0,  we  can  integrate  the  differentiated  form  obtained  by  combining  the  identity 
dS=C^dT/T+(dp/dT)vdV  and  Eq.  (2.3 . 1 1 7),  to  get 

TV”  =  expJdS/C^  .  (2.3.118) 


It  follows  from  Eq.  (2.3.1 16)  that  we  caimot  perform  the  integration  in  Eq.  (2.3.1 18)  until  is 
know  as  a  function  of  S.  We  simply  bypass  this  difficulty  here  by  assuming  that  is  a  constant. 
Then  E,  satisfies  the  equation  dE^ldS^e^'^"  which  on  integration  yields  the  equation 
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E^(S)  =  C^TV"+E°  .  (2.3.119) 

Combining  Eqs.  (2.3.1 13)  and  (2.3.1 19)  and  setting  n=Y-l  gives  the  p=p(T,V)  EOS  as 

pV=(y-l)C^T  +  E°/V''-^  (2.3.120) 

with  El  a  parameter  to  be  determined. 

We  finally  assume  that  our  polytropic  fluid  has  a  mass  m,  set  V=mv,  and  rewrite  Eq. 
(2.3.120)  in  terms  of  specific  quantities  as 

A 

pv  =  (y-l)cT  +  i(-)T^'*  (2.3.120)' 

V 

with  c^=Cv/m  and  a  constant  A  =  (pv-(y-\)c^f) .  Equation  (2.3.120)'  can  be  calibrated  for 
detonation  products  by  using  known  values  of  p  ,  v  ,  f  ,  and  c^at  the  CJ  point  to  evaluate  the 
constant  A . 

2.4.  Steady-State,  Nonideal  Detonation  in  an  Explosive 

2.4.1.  Composition  Containing  A1 

Our  thermodynamic  models  for  explosive  compositions  presented  in  Section  2.1  and  2.2 
are  both  based  on  the  tenet  that  a  satisfactory  treatment  of  nonideal  detonation  in  an  explosive 
composition  containing  aluminum  (Al)  must  take  accoimt  of  the  reaction  of  the  explosive  grains 
and  Al  particles  in  the  reaction  zone  (RZ)  neighboring  the  shock  front.  This  tenet  leads  to  two 
limiting  cases  of  nonideal  detonations,  that  propagate  at  constant  velocity  and  consequently  are 
supported  by  steady-state,  reaction  zones  (SSRZ).  In  one  of  these  cases,  the  time  scale  for  the 
reaction  of  the  Al  with  the  explosive’s  products  is  much  larger  than  the  time  scale  for  the 
decomposition  of  the  explosive  grains:  all  the  Al  passes  through  the  SSRZ  before  it  reacts  with 
the  product  gases  in  the  ensuing  release  wave.  In  the  other  case,  the  time  scale  for  the  reaction  of 
Al  with  the  explosive’s  products  is  significantly  smaller  than  the  time  scale  for  the 
decomposition  of  the  explosive  grains:  all  the  Al  reacts  in  the  SSRZ;  the  release  wave  can  be 
modeled  as  an  isentrope.  Consequently,  a  constant- velocity,  nonideal  detonation  will  in  general 
be  represented  by  an  intermediate  case  with  Al  reacting  in  both  the  SSRZ  and  the  release  wave, 
and  rate  functions  for  the  decomposition  of  the  explosive  and  the  reactions  of  Al  must  be  known 
before  a  thermo-hydrodynamic  model  for  such  a  detonation  can  be  used  to  calculate  its 
properties.  For  these  reasons,  we  performed  two  tasks  to  provide  a  more  definitive  backgroimd 
for  nonideal  detonation  and  thereby  develop  a  better  appreciation  of  its  complications.  In  the 
first  task,  we  developed  a  prototype  model  for  SS,  nonideal  detonation,  and  in  the  second,  a 
prototype  model  for  burning  explosive  grains.  The  work  performed  on  the  first  task  will  be 
presented  in  this  section  but  the  work  performed  on  the  second  task  will  be  presented  in 
Appendix  A. 

To  provide  background  for  our  treatment  of  SS,  nonideal  detonation,  we  first  present  the 
equations  governing  a  SS  detonation  wave  and  then  consider  the  relationship  between  its  RZ 
length  and  corresponding  reaction  time. 


36 


ENIG  TR  95-1 


2.4.2.  Equations  for  a  Steady-State  Reaction  Zone 

Here,  we  will  use  the  subscripts  x  and  p  to  denote  the  explosive  and  its  products,  the 
subscript  H  to  denote  the  shocked  state,  and  either  a  superscript  or  a  subscript  0  to  denote  the 
standard  state  of  the  explosive  and  the  detonation  products,  i.e.,  at  pressure  po  and  temperature 
Tq.  Recalling  that  the  Rankine-Hugoniot  (RH)  jump  conditions  are  first  integrals  of  the 
differential  equations  expressing  the  conservation  laws  when  the  flow  has  attained  a  SS,  we  write 
the  equation  governing  our  SSRZ  as 


p(D-u)=poD, 

(2.4.1) 

P=PoD^(l-v/Vo), 

(2.4.2) 

e-ex=pvo(l-v/vo)/2, 

(2.4.3) 

with  D  and  u  used,  respectively,  to  denote  the  detonation  velocity  and  the  particle  velocity,  and 
Po=0.  Eq.  (2.4.1)  expresses  the  conservation  of  mass,  Eq.  (2.4.2)  expresses  the  conservation  of 
mass  and  momentum,  Eq.  (2.4.3)  expresses  the  conservation  of  mass,  momentum,  and  energy, 
and  the  equations  governing  the  shock  at  the  front  of  the  SSRZ  are  obtained  by  putting  a 
subscript  H  on  the  variables  p,u,p,v  and  e  in  Equations  (2.4.1),  (2.4.2),  and  (2.4.3). 

To  continue  this  treatment  of  our  SSRZ,  we  introduce  the  sound  speed  c,  define  the 
adiabatic  index  k  by  the  equation  c^=kpv,  and  recall  that  the  SSRZ  terminates  at  the  Chapman- 
Jouquet  (CJ)  state  where  the  flow  become  sonic.  Using  the  subscript)  to  denote  the  CJ  state,  this 
sonic  condition  is  expressed  by  the  equations, 

D-Uj=Cj,  (2.4.4) 

The  combination  of  the  equations  pjCj=poD,  c?=kjPjVj,  and  Eq.  (2.4.2)  written  for  the  CJ  state,  then 
gives  the  well  known  equations  for  v,u,c,  and  p  in  the  CJ  state: 

Vj  /vo=kj/(kj+l),  (2.4.5) 

kjUj=Cj=Dkj/(kj+l),  (2.4.6) 

Pj=PoD7(kj+l).  (2.4.7) 

We  finally  note  that  the  relationship  between  D  and  the  specific  chemical  energy  q 
liberated  by  a  particle  passing  through  the  SSRZ  is  derived  by  combining  the  e=e(p,v)  equation 
of  state  for  the  detonation  products  with  the  following  equation: 

ej-e®=DV2(kj+l)2  (2.4.8) 


obtained  by  combining  Eqs.  (2.4.5)  and  (2.4.7)  with  Eq.  (2.4.3)  written  for  the  CJ  state.  Having 
derived  the  RH  jump  conditions  for  SS  flows,  it  is  convenient  to  consider  how  RZ  length  and 
reaction  time  are  related  in  a  SS  detonation. 
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2.4.3.  The  Relation  Between  Reaction  Zone  Length  and  Reaction  Time  in  a  Steady 
State  Detonation  Wave 

Our  purpose  here  is  to  derive  equations  that  place  upper  and  lower  limits  on  the  SSRZ 
length  L  when  the  detonation  velocity  D  is  known  and  the  reaction  time  t^  has  been  estimated. 
We  accordingly  consider  a  one-dimensional,  SS,  ZND  (Zeldovich,  von  Neumann,  Doering) 
detonation  wave.  With  our  present  notation,  the  equation  for  the  RZ  length  can  be  written  as 

L  =  Dt^  -  fu  dt  ,  (2.4.9) 

0 


where  t  denotes  the  time. 

Because  u  drops  linearly  Avith  p  along  the  Rayleigh  line  representing  the  SSRZ,  realistic 
equations  for  lower  and  upper  limits  of  L  can  be  obtained  by  setting,  respectively,  u=Uh  and  u=Uj 
in  Eq.  (2.4.9).  To  derive  the  equation  for  the  lower  limit  L.,  we  assume  that  shock  velocity  U,Uh 
Hugoniot  curve  for  the  explosive  is  linear,  U=a+Bu,  and  use  the  value  of  UH=(D-a)/B  in  Eq. 
(2.4.9)  to  obtain  the  equation 


L  =D(l-(D-a)/DB)tR.  (2.4.10) 

To  derive  the  equation  for  the  upper  limit  L+,  we  recall  that  Uj=D/(kj+l),  where  k  denotes  the 
adiabatic  index  (dlnp/dlnv),,,  and  use  this  value  for  Uj  to  obtain  the  equation 

L^  =  kjDtR/(kj+l)  (2.4.11) 

with  the  values  of  a=2.90  mm/psec  and  B=1 .59  for  AP  given  by  SEP"  and  the  values  of  D=6.6 
mm/psec  and  kj=3.45  for  AP  obtained  from  a  TIGER  code  calculation.  Equations  (2.4.10)  and 
(2.4.1 1)  for  AP  become 


L=4.29tR  (2.4.10)' 

L^  =  5.16tR.  (2.4.11)' 

Thus  if  tR=0.5  psec,  the  upper  and  lower  limits  for  the  length  of  the  RZ  in  AP  are,  respectively, 
2.58  mm  and  2.14  mm. 

It  is  now  convenient  to  formulate  our  prototype  model  for  SS  nonideal  detonation. 

2.4.4.  A  Prototype  Model  for  Nonideal  Steady  State  Detonation 

For  convenience  in  formulating  a  prototype  model  for  nonideal  SS  detonation,  we  again 
consider  the  C^HbNPdAlf  composition  discussed  in  Section  2.2.  We  adopt  the  reaction  scheme 
formulated  for  this  composition,  but  for  the  sake  of  tractability,  make  more  simplistic 
assumptions  about  the  condensed  components.  We  formulate  the  simplest  treatment  of 
detonation  supported  by  (R.3)  and  (R.7)  by  assuming  the  following: 
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A.7  The  explosive  X=CaHbNgOd  is  a  polytropic  explosive  with  the  same  index  as  the 
polytropic  products. 

A.8  Both  A1  and  AI2O3  can  be  treated  as  incompressible  solids,  and  the  reaction  between 
A1  and  HjO  produces  only  one  form  of  AI2O3.  With  these  assixmptions,  the  e-relationship 
we  previously  developed  for  our  CaH^NeOdAlf  composition  takes  a  much  simpler  form, 
and  the  subscript  and  superscript  used  to  denote  different  forms  of  AI2O3  can  be 
eliminated. 

To  be  more  definitive,  we  denote  the  specific  volumes  of  X,  Al  and  AI2O3  by  v,^,  Vj ,  and 
V2  and  make  two  observations.  The  first  observation  is  that  the  A^^  and  A°  terms  disappear  from 
Eq.  (2.2.23)  and  (2.2.25)  because  from  A.7.  The  second  is  that  Eq.  (2.2.24)  becomes 

Eg=0  because  gx(V)=0  from  A.7  and  gAi(y)=gi(V)=0  from  A.8.  In  this  case,  Eq.  (2.2.25)  can  be 
rewritten  in  a  simpler  form  as 


y-l  y-l 


(2.4.12) 


by  setting  el=e^  ,  ,  q^^q^^  ,  and  l:s.B=B^-B^ 


It  is  now  convenient  to  derive  the  relationship  between  the  adiabatic  index  k  and  y.  We 
consequently  use  the  thermodynamic  identiy  for  the  sound  speed  and  c^=kpv,  which  yield 


V  de 
dp 


(2.4.13) 


and  equations  for  the  partial  derivatives  of  e  derived  from  Eq.  (2.4.12)  to  obtain  the  equation 


k-- 


yv 


(2.4.14) 


We  now  use  Eqs.  (2.4.12),  (2.4.14),  and  (the  equations  expressing  the  CJ  conditions)  Eqs. 
(2.4.5)-(2.4.8)  to  derive  equations  for  steady  state  detonation  parameters  and,  thereby, 
demonstrate  the  major  problem  we  encounter  in  modeling  nonideal  detonations. 


2.4.5.  The  Equations  for  Steady  State  (SS)  Detonation  Parameters 


To  derive  equations  for  SS  detonation  parameters,  we  first  combine  Eqs.  (2.4.14)  and 
(2.4.5)  to  obtain  an  equation  for  the  CJ  volume  Vj.  Because  Eq.  (2.4.14)  contains  the  reaction 
coordinate  we  are  immediately  faced  with  the  main  problem  for  this  type  of  nonideal 
detonation,  namely,  that  of  determining  how  far  the  reaction  of  Al  with  the  detonation  products 
has  progressed  at  the  CJ  point.  To  account  for  but  not  solve  this  problem,  we  assume  that  ^2j  is  a 
parameter  to  be  determined  and  write  our  equation  for  Vj  from  Eq.  (2.4.14)  and  (2.4.5)  as 


yvo-(^i+^^^2,) 

y+1 


(2.4.15) 
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The  equation  for  (kj+l)'*  then  follows  by  combining  Eqs.  (2.4.5)  and  (2.4.15)  as 

1  ^0  ’^0 


(2.4.16) 


k.+\ 

j 


Y+1 


and  Eqs.  (2.4.16)  and  (2.4.7)  give  the  equation  for  the  CJ  pressme  as 

'>  Lb  y- 


Pj  = 


^0  ^0 


(2.4.17) 


Y+1 


Equations  (2.4.8),  (2.4.12),  (2.4.15),  (2.4.16),  and  (2.4.17)  then  give  the  equation  for  the 
detonation  velocity  as  follows.  We  first  combine  Eqs.  (2.4.8)  and  (2.4.16)  to  obtain  the  equation 

Lb,  ^ 


D‘ 


(2.4.18) 

2(Y+1)^  Vq  V, 

We  then  combine  Eq.  (2.4.12),  with  ^,=1  and  Eqs.  (2.4.15)  and  (2.4.17)  to  obtain  the 

equation 

[1+^+My' 

_  s  .  Vo  Vo 

(Y*l)'  Y-l 

The  combination  of  Eqs.  (2.4.18)  and  (2.4.19)  then  gives  the  equation  for  D  as 


(2.4.19) 


= 


2(Y'-i)(gi+g2y 

V  V  ^ 


(2.4.20) 


Eqs.  (2.4.15)-(2.4.17)  and  Eq.  (2.4.20)  show  explicitly  that  the  value  of  must  be  known  before 
SS  detonation  parameters  can  be  calculated  when  all  the  other  properties  of  the  explosive 
composition  are  known.  They  also  reduce  to  the  equations  for  the  two  limiting  cases  of  nonideal 
detonation  discussed  at  the  begiiming  of  this  section  as  follows.  When  the  reaction  of  the  A1 
with  the  explosive’s  products  is  very  slow  compared  to  the  decomposition  of  the  explosive,  we 
can  obtain  the  equations  for  the  SS  detonation  parameters  by  setting  m  the  above 
equations.  When  the  reaction  of  the  A1  with  explosive’s  products  is  appreciably  faster  than  the 
decomposition  of  the  explosive,  we  can  obtain  the  equations  for  the  SS  detonation  parameters  by 
setting  ^2j=l  in  the  above  equations. 


In  general,  it  is  clear  that  an  expression  for  the  rate  of  reaction  of  the  A1  with  the 
detonation  products,  d^2/<it,  must  be  known  before  the  value  of  ^2j  can  be  calculated  and  used  in 
Eqs.  (2.4.15)-(2.4.17)  and  (2.4.20)  to  calculate  the  SS  detonation  parameters  for  a  nonideal 
detonation. 
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3.  THE  AP/AI/EXPLOSIVE  GAS  REACTIONS 

3.1  The  A1  Reaction  Model 

The  reaction  of  aluminum  (Al)  with  air  imder  normal  ambient  conditions  is  a  well  known 
phenomenon.  Under  these  conditions  a  surface  of  pure  Al  will  become  coated  with  a-Al203 
(corundum),  which  will  grow  in  thickness  to  a  few  tens  of  A  in  a  few  days.  This  process  is 
electrolytic  in  nature  has  been  extensively  studied  and  is  fairly  well  understood.^^'^® 

In  this  process,  a  battery  is  formed  with  a-Al203  playing  the  role  of  solid  electrolyte.  The 
AI/AI2O3  interface  is  the  battery  anode,  and  the  air  (02)/Al203  interface,  the  cathode.  At  the 
anode,  AU^  ions  (cations)  are  absorbed  into  AI2O3  and  migrate  toward  the  cathode.  Likewise,  at 
the  cathode,  O"  ions  (anions)  may  be  absorbed  and  migrate  toward  the  anode.  AI2O3  is  created 
through  the  interactions  of  the  cations  and  anions. 

One  can  view  the  current  inside  the  AI2O3  in  terms  of  the  movement  of  cation  vacancies 
(hereafter  called  vacancies)  and  vacancies  in  the  electron  band  of  the  oxide^*  (hereafter  called 
holes).  Let  p+  and  p.  be,  respectively,  the  chemical  potentials  of  the  cation  vacancies  and 
electron  holes,  and  let  the  AI2O3  be  an  infinite  slab  bounded  by  two  parallel  planes  a  distance  I 
apart.  Let  x  be  the  perpendicular  distance  from  one  plane  and  x=0,  define  the  the  anode  surface 
(the  AI/AI2O3  interface).  The  cathode  surface  (the  AI2O3/O2  interface)  is  then  at  x=/. 

In  the  steady  state  of  current  flow,  the  number  of  cation  vacancies  passing  through  a  plane 
perpendicular  to  x  must  be  balanced  by  three  times  the  number  of  holes  passing  in  the  opposite 
direction. 


The  system  is  functionaly  a  galvanic  cell  and  obeys  the  same  electrochemical  laws.  One 
therefore  expects  charge  neutrality  (or  near  neutrality)  throughout  the  body  of  the  electrolyte  and, 
in  addition,  that  the  gradient  of  the  chemical  potential  be  balanced  by  the  electrostatic  force^^’^^ 
that  is. 


^^(P+  +  3p.) 
dx 


3ee  , 


(3.1.1) 


where  e  is  the  electron  charge  and  e  the  electric  field.  The  current  density)  in  the  cell  is  the  sum 
of  two  terms,  j+  and  j.,  which  are  equal.  The  current  density  j+  is  that  for  cation  vacancies,  and  j. 
for  anion  vacancies.  The  field  e  is  related  to  j  via  Ohm’s  law,  i.e., 

JP  =  e.  (3.1.2) 


In  Eq.  (3.1.2),  the  resistivity  p  is  the  sum  of  two  resistivities,  p+  and  p.,  one  for  cations  (or 
vacancies)  and  one  for  electrons  (or  holes).  These  have  eorresponding  conductivities  0+  and  a  ,so 
Eq.  (3 . 1 .2)  can  be  written 


0^+0 

e  =J - 

0,0. 


(3.1.3) 


Using  this  to  substitute  for  €  in  (3.1.1),  one  can  solve  for  current  density  and  find 
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j  = - • 

3e(a^+o_)  dx 

The  growth  of  electrolyte  (the  creation  of  AI2O3)  is  simply  related  to  j  by 

dl  _  j^e 
dt  6e 


(3.1.4) 


(3.1.5) 


in  which  t  is  time,  j/6e  is  half  the  number  of  cations  plus  corresponding  holes  passing  through  a 
unit  area  per  unit  time  and  Qe  is  the  volume/molecule  of  AI2O3.  Note  that  for  a+  and  o.  one  may 
use  the  basic  semiconductor  relationships 

(J.  =  3en,il;,  (3.1.6) 


o-  =  enjf_  ,  (3.1.7) 

in  which  n±  and  1);^  are,  respectively,  the  number  densities  and  mobilities  for  vacancies  and  holes. 
Because  of  the  neutrality  requirement,  n+  =  n ,  and  as  will  be  shown  further  on,  the  derivative  in 
Eq.  (3.1.4)  should  equal  (4kT/iv)dn+/dx,  where  k  is  the  Boltzmarm  constant  and  T  is  temperature 
in  °K.  With  these  considerations,  using  (3.1 .6)  and  (3.1 .7)  for  0+  and  o.,  (3 . 1 .4)  becomes 

4A:rv|f^l(;  dn^ 

J  =  •  (3-1-8) 


One  may  properly  assume  ilt.»tlr+  and  integrate  both  sides  of  (3.1.8)  between  x=0  and 
x=l  to  obtain 


j  = 


«  m-«iO) 


(3.1.9) 


This  expression  for  j  contains  calculable  quantities  on  the  right  side  and  may  be  substituted  into 
Eq.  (3.1.5)  to  obtain  the  rate  equation 


dt  I 


where  the  rate  constant  K  is 


is:  =  4Q 


«.(/)-«  .(0) 


Equation  (3.1.10),  when  integrated,  gives  the  parabolic  rate  law 

=  2Kt  . 

One  should  note  that  in  terms  of  the  ionic  conductivity  o+,  K  has  the  form 

oTO-oYO) 

K  ^  AQkT  ^  ^  . 


(3.1.10) 


(3.1.11) 


(3.1.12) 


(3.1.13) 


These  last  four  expressions  are  the  basic  model  for  the  reaction  rate  of  Al,  which  will  be 
discussed  and  developed  further  in  succeeding  sections. 
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3.2  A  Plausibility  Argument  Based  on  Existing  Data 

Experimentally,  A1  has  been  reacted  with  oxygen  both  in  explosive  detonations^*"''^  and  in 
controlled  laboratory  experiments.^  Experiments^^  with  flat  slabs  of  pure  A1  in  gaseous  O2  (at  76 
torr)  have  been  performed  to  temperatures  above  800°K.  The  reaction  product  is  on  admixture  of 
crystalline  and  amorphous  Y-AI2O3,  and  the  reaction  rates  even  at  the  elevated  temperatures  are 
quite  slow.  The  thickening  of  the  products  typically  occurs  at  rates  of  a  few  microns  per  hour. 
The  reaction  rate  constant  grows  with  T,  increasing  by  a  factor  of  about  8  between  750“K  and 
850°K. 


Growdh  rates  follow  a  parabolic  rate  law  imtil  some  limiting  film  thickness  is  reached. 
These  thicknesses,  typically  150A  to  200A,  increase  linearly  with  temperature  at  very  roughly 
0.5A  per^K. 

There  are  a  number  of  relevant  experiments  involving  aluminized  explosives.  These 
produce  reaction  results  varying  from  partial  to  more  or  less  complete  reaction  of  Al.  Some  of 
this  work  involve  AP  and  Al  exclusively,^*’®’''^  and  in  one  effort  a  number  of  test  samples  are 
close  to  the  theoretical  maximum  density  (TMD).''^  In  a  few  of  the  latter  experiments,  complete 
or  nearly  complete  conversion  of  Al  to  AI2O3  is  observed.  Since  the  mean  diameter  of  Al 
particles  lies  between  25p  and  30p  and  the  experimental  time  scale  is  of  the  order  of  lOpsec,  the 
implied  linear  burning  rate  of  Al  is  in  the  range  of  hundreds  of  cm  per  second.  This  is  nine 
orders  of  magnitude  greater  than  the  p/hour  rate  discussed  above  in  the  oxidation  experiments 
and  implies  a  rate  constant  1 8  orders  of  magnitude  greater  than  those  observed.  In  the  oxidation 
experiments,  the  reaction  rate  with  temperature  is  seen  to  increase  by  a  factor  of  about  8  over 
lOO^K.  At  this  rate  of  eight-fold  increase,  a  factor  of  10'*  is  achieved  after  a  temperature  increase 
of  2000°K.  Since  the  maximum  experimental  oxidation  rates  for  pure  Al  are  observed  at  850°K, 
the  rates  observed  in  explosive  reactions  should  be  reached  at  about  2850°K,  which  is  in  the 
proper  range  for  explosive  reaction  products.  In  other  words,  a  complex  temperature 
extrapolation  of  the  oxidation  data  moves  it  into  the  explosive  regime,  and  this  without  regard  to 
the  O2  pressure,  which  in  an  explosive  is  10®  times  greater  than  the  76  torr  in  the  oxidation 
experiments,  further  increasing  the  reaction  rate.  This,  however,  cannot  be  said  of  the  limiting 
thicknesses  of  AI2O3  observed  in  these  oxidation  experiments.  Linear  extrapolation  of  these  data 
to  2850“K  indicates  an  increased  thickness  of  about  450A,  which  amounts  to  only  a  few 
hundredths  of  1  p.  Since  explosive  experiments  do  not  observe  such  small  limits,  the  thickness 
limits  are  either  drastically  affected  by  pressure  or  very  nonlinear  in  temperature  or  both. 

Another  set  of  seemingly  relevant^®  data  involves  the  P-aluminas,  which  exhibit  the 
phenomenon  of  anomalously  high  ionic  conduction.  In  structure,  these  are  basically  slabs  of 
crystalline  y-aluminum  (about  lOA  thick)  separated  by  a  layer  or  layers  of  monovalent  metal 
oxide.  Within  this  layer  are  “tunnels”  through  which  ions  can  move  rather  fireely.  These 
materials  exhibit  conductivties  as  high  as  0.1  (Qcm)"'  at  1000°K.  Using  this  value  in  (3.1.13), 
taking  Qg~4.3xl0"^®,  and  assuming  a+(0)=0,  one  can  compute  K  and  substitute  into  Eq.  (3.1.12) 
to  find  /  at  10  psec.  This  yields  /«70p,  a  value  also  consistent  with  the  results  of  explosive 
experiments. 

All  the  above  data  can  be  utilized  to  piece  together  some  plausible  explanation  of  events 
in  Al/explosive  reactions.  First  of  all,  x-ray  analysis  of  explosive  products  reveals''®  a-Al203  and 
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not  the  Y-AI2O3  observed  in  low  pressure  oxidation  experiments  (T<850°K).  This  makes  some 
sense  because  Y-AI2O3  is  known  to  transform  to  a-Al203  at  elevated  temperature,  and  moreover, 
a-Al203  is  the  denser  of  the  two  crystalline  phases  and  should  be  preferred  at  very  high  pressures. 
The  x-ray  analyses,  however,  are  qualitative  and  do  not  exclude  amorphous  AI2O3  in  the 
explosive  products,  nor  do  they  exclude  some  P-like  form  of  Y-AI2O3,  which  may  accrue  due  to 
the  presence  of  monovalent  metal  impurities.  The  current  model  will  assume,  therefore,  that 
both  a-Al203  and  some  amorphous  form  of  Y-AI2O3  are  produced  in  the  explosive  gases.  The 
general  AI2O3  products  will  be  considered  to  contain  tunnel  structures  through  which  AT^  cations 
can  readily  pass  and  from  which  O”  are  excluded  because  of  their  much  larger  ionic  radius  (2.5 
times  that  of  AE^).  The  very  last  assumption  is  made  with  the  knowledge  that  if  significant 
reaction  occured  at  the  Al/oxide  interface  it  would  have  little  effect  on  the  explosive  gases 
because  of  the  time  necessary  for  the  heat  of  reaction  to  conduct  outward  through  the  AI2O3 
product  layer.  This  time  is  about  where  P  is  reciprocal  diffusivity.  For  AI2O3,  P«900 
sec/cm^,  so  that  the  time  for  the  reaction  heat  to  pass  through  1  p  of  product  film  is  9  psec  or 
about  the  full  length  of  an  explosive  experiment.  Eventually,  the  energy  escapes,  but  there  would 
be  little  effect  on  the  overall  release  isentrope  of  the  explosive  gases  and,  therefore,  on 
preformance. 

The  only  term  in  Eq.  (3.1.11)  that  can  account  for  the  huge  increase  in  K  under  explosive 
conditions  is  the  cationic  mobility  i|;+.  In  Y-AI2O3  (and  presumably  in  a-Al203),  this  has  an 
Arrhenius-like  temperature  dependence,  which  cannot  explain  the  increase  in  question  unless  the 
activation  energy  itself  decreases  drastically.  Since  the  activation  energy  is  a  measure  of  the 
energy  barrier  a  cation  must  overcome  to  change  sites,  one  would  expect  an  increase  in  this 
barrier  energy  with  pressure  (decreasing  volume)  as  opposed  to  a  decrease.  If  barrier  height, 
however,  did  decrease  in  some  unusual  way  with  decreasing  volume,  the  relative 
incompressibility  of  a-Al203 ,  which  is  gem-like,  would  prevent  much  change  from  occuring. 

The  postulation  of  tunnels  as  cation  conduits  seems  necessary  to  obtain  lower  activation  energies 
for  cation  motion  because  of  the  relative  incompressibility  of  a-Al203,  and  because  the 
temperature  extrapolation  discussed  above  is  highly  nonlinear  (and  not  very  credible)  and  leads 
us  to  the  very  highest  part  of  the  range  of  detonation  temperatures.  A1  is  known  to  react 
agressively  with  oxygen  at  much  lower  temperatures  (but  above  about  1000“K). 

In  siunmary,  interpretation  of  the  data  discussed  above  leads  to  the  following  preliminary 
reaction  model: 

1 .  Tunnel  structures,  which  conduct  AF^  cations,  are  assumed  to  form  naturally  with  the 
AI2O3  reaction  products  in  detonated  Al/explosive  compositions. 

2.  Anions  are  excluded  from  the  conducting  tunnels 

3.  The  solid  reaction  products  are  assumed  to  be  an  admixture  of  a-Al203,  amorphous 
AI2O3,  and,  possibly,  some  P  form. 

3.3  Equilibrium  at  the  Cathode  and  Anode 

At  the  anode,  A1  atoms  enter  AI2O3  as  AF^  ions  giving  3  electrons  to  the  AI2O3  bands, 
i.e.,  within  the  AI2O3, 
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Al  ^  Ar^  +  3e  .  (3.3.1) 

At  the  cathode,  reacts  with  Oj  to  form  AljOj,  i.e., 

6e+2Al^^  +  ~^2  ^  Al^O^  +  2F^  +  6V_  ,  (3.3.2) 

where  V+  represents  a  cation  vacancy  and  V.  a  hole  in  the  electron  bands.  The  equilibrium  in 

these  reactions  may  be  used  to  find  the  quantity  n+  in  Eq.  (3. 1 .1 1).  To  simplify  the  process, 
subscripts  will  be  assigned  to  connect  the  properties  of  a  constituent  with  the  constituent.  The 
subscript  e  will  detonate  properties  of  AI2O3  (such  as  Q^,).  In  addition,  a,  x,  c,  and  q  will 
designate  properties  of  Al,  Oj,  Af^  and  electrons,  respectively.  As  before,  +  and  -  will  describe 
properties  of  vacancies  and  holes. 

At  equilibrium,  for  the  reaction  Eq.  (3.3.1),  one  can  write^’’'*'* 

,  (33.3) 

where  the  p  are  the  indicated  chemical  potentials.  Likewise  for  reaction  Eq.  (3.3.2),  one  has 

>6n.  .  (3.3.4) 

From  statistical  mechanics,^’  ''''  one  knows  that  for  a  unit  of  the  constituent  k, 

H  ’  (3.3.5) 


where  q^  is  the  partition  function  per  unit,  is  the  number  of  the  units  that  are  indistinguishable 
from  one  another,  andN^^l  if  they  are  all  distinguishable. 


In  this  regard,  both  the  atoms  of  Al  in  the  metal  and  the  molecules  of  AI2O3  in  the 
electrolyte  are  completely  distinguishable  because  of  their  fixed  positions.  This  also  applies  to 
the  electrons  in  the  electrolyte  energy  bands,  which  are  distinguishable  by  their  quantum 
numbers.  Therefore, 


N  ^  N  =  N  =  1 

a  ‘  c  q 


(3.3.6) 


On  the  other  hand,  vacancies,  holes,  and  O2  molecules  are  completely  indistiguishable 
because  they  have  neither  fixed  positions  nor  quantum  numbers.  Therefore,  N,j,  N+  and  N  are 
simply  equal  to  the  respective  numbers  of  these  constituents. 


The  situation  for  cations  travelling  in  tunnels  is  more  complex  because  cations  within  a 
given  tunnel  are  indistinguishable,  but  cations  in  different  tunnels  can  be  distinguished.  This 
subject  is  treated  next. 

3.4  The  Partition  Function  Per  Cation 


In  travelling  from  anode  to  cathode  a  cation  will  be  assumed  to  traverse  a  tunnel  of 
(mean)  length  /j.  Since  the  system  is  expected  to  be  at  high  temperature,  it  is  described  by  the 
semiclassical  grand  partition  function'’'* 
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Qr 


—  I  ..  I  e  ...  dpjq^ 


(3.4.1) 


h^\w,\r 

Here  h  is  Planck’s  constant,  p,q  are  the  momentum  and  space  coordinates,  respectively, 
P=l/kT,  w  is  the  total  number  of  cations,  w,  is  the  number  of  cations  per  tunnel,  Nj  is  the  number 
of  tunnels,  and  H  is  the  Hamiltonion.  How  Eq.  (3.4.1)  is  evaluated  depends  on  the  configuration 
of  the  tunnels.  If  they  are  very  wide  and  high  compared  to  the  Al^^  diameter  (~  1  A),  then  the 
cation  can  be  viewed  as  bouncing  from  the  tunnel  walls  while  traversing  the  length.  If  the  width 
and/or  height  are  relatively  narrow,  then  the  cation  may  be  vibrating  transversely  in  one  or  both 
directions  while  travelling  the  length.  Along  the  tunnel  length,  there  may  be  barriers  to  motion 
that  must  be  overcome.  Let  the  distance  between  barriers  be  s  and  the  distance  across  a  barrier 
be  Sfc.  Also,  let  the  barrier  height  be  e^.  Assume,  as  well,  that  the  tunnel  is  narrow  in  either  or 
both  of  the  two  transverse  directions  and  that  the  ion  oscillates  harmonically  in  these  directions 
with  frequencies  and  for  the  height  and  width,  respectively.  The  number  of  cations  Wj 
within  a  tunnel  is  expected  to  be  small  relative  to  the  number  of  sites  of  length  s+s,,.  This  being 
the  case,  the  Hamiltonion  H^  of  a  single  cation  can  be  written  as 

,2  „ 


H 


_  P 


2m 


4lt^V^— 

9  T 


+  e„ 


(3.4.2) 


where  m^  is  the  cation  mass,  e®  is  the  potential  of  the  cation  within  the  tunnel  when  in  the  region 
of  length  s  between  barriers,  the  function  g(q)  is  zero  within  the  region  s  and  equals  one  in  the 
barrier  region  s,,,  and  the  oscillation  frequecy  is  assumed  to  be  zero  if  q  is  in  a  direction  of 
large  tunnel  width  or  height. 


Using  Eq.  (3.4.2)  in  Eq.  (3.4.1)  for  a  single  cation,  the  integration  can  be  performed  over 
the  three  space  and  three  momentum  coordinates,  yielding  the  integral 


/  -(iTim  kT) 


fkT  ,kT  i-f'h 

( — )  (— )  X. 


(3.4.3) 


In  this  expression,  the  indices  V^,  are  unity  if  the  height  h  or  width  w  of  the  tunnel  are  small 
enough  to  cause  transverse  harmonic  vibrations  of  the  cation,  and  are  zero  otherwise.  One 
should  also  note  that  the  ratio  //(s,+Si,)  is  just  the  number  of  sites  of  length  s+s,,  within  a  tunnel. 
This  ratio  carries  the  integration  over  a  single  site  through  the  entire  tunnel  length. 


The  integral  Eq.  (3.4.3)  is  over  the  coordinates  of  a  single  cation.  If  integration  is 
performed  for  all  cations  within  a  tunnel  one  obtains  (Z^)"^'.  Doing  this  for  all  tunnels,  one 
obtains  (/^)^''^'for  the  multiple  integral  in  Eq.  (3.4.1).  Since 

^  (3.4.4) 

then  Eq.  (3.3.7)  ean  be  rewritten  as 


Qc 


(3.4.5) 


Using  the  Stirling  approximation 
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"W  I - 

w^\  ~  e 


(3.4.6) 


and  remembering  that  N,Wt=w,  one  finds 


.N,  w 


(w^!)  '  ~  Wf  e  (2tiw^ 

The  Hemholtz  energy  for  the  cations  A^,  is  given  by 

-  -kT  \nQ^  . 

Substituting  Eq.  (3.4.7)  into  Eq.  (3.4.5)  for  Q^,  one  finds  for  A^,, 

=  -A:rvi'[ln -  +  1- — -  ln(2Ti:vr^)]  . 

h^w  2w 


(3.4.7) 


(3.4.8) 


(3.4.9) 


Since  N/w«l,  one  can  neglect  the  last  two  terms  in  the  brackets,  and  A,,  becomes 

A  =  -kTw  In -  . 

h 

Using  Eq.  (3.4.10),  one  can  obtain  the  chemical  potential  for  cations  from  the  relation 


(3.4.10) 


which  on  neglecting  small  terms  yields 


p  =  -kT  In- 


(3.4.11) 


(3.4.12) 


Relating  Eq.  (3.4.12)  to  Eq.  (3.3.5),  one  can  see  that 


(3.4.13) 


•  (3.4.14) 

t 

This  completes  the  statistical  description  of  cations;  the  development  is  necessary  because 
it  is  specific  to  the  model  and  has  not  previously  been  done.  The  partition  functions  of  the 
remaining  constituents  are  better  understood  and  are  given  below. 


3.5.  The  Evaluation  of  the  Remaining  Partition  Functions:  q^,  q^,  q^,  q+,  q.,  q^ 

Subsequent  to  shocking,  the  unreacted  A1  is  expected  to  be  above  its  Debye  temperature. 
So  for  q^,  we  find 


4 

9a  tg  ; 


(3.5.1) 


where  0^  is  the  Debye  temperature  divided  by  the  cube  root  of  the  base  of  the  natural  logrithm 
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and  6°  is  the  potential  energy/atom  of  A1  in  the  metal  at  0°K  (a  function  of  volume  only)  plus  the 
energy  required  to  transfer  three  electrons  to  the  AljOj  band.  For  AI2O3,  a  single  Debye 
temperature  will  be  assigned  to  the  vibrational  modes  per  degree  of  freedom  of  the  molecule  (3 
associated  with  translation  and  3  with  rotation).  For  these  acoustic  modes,  T  is  again  expected  to 
be  greater  than  the  assigned  Debye  temperature.  For  each  of  the  nine  remaining  vibronic  modes, 
an  Einstein  temperature  Og  is  assigned.  This  gives  for  qg. 


T  ^ 


2T 


-]  ^ 


-P[e‘+e!] 


1  -e 


(3.5.2) 


where  6g  is  defined  in  the  same  manner  as  for  Al,  the  Oj  are  the  nine  Einstein  temperatures,  is 
the  binding  energy  of  the  AI2O3  molecule,  and  e®  is  the  potential  of  the  molecule  in  the  AI2O3 
lattice. 


The  electron  partition  function  q  is  quite  simply 


-Pe" 

9,  =  e  ’ 


(3.5.3) 


in  which  6°  is  the  energy  binding  the  electron  in  the  AI2O3  band.  The  functions  q+  and  q_  are 
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9.  =  Ne 


-Pe! 


(3.5.4) 


and 


9-  =  Ne 


-Pe“ 


(3.5.5) 


The  quantities  e+  and  e?  are  the  negatives  of  e®  and  e®,  respectively. 
The  partition  function  q,^  for  O2  gas  is 


9,  =  ( 


iTZm  kT  - 

Ve 


Orb  0., 

^--^2 


(3.5.6) 


where  m^  is  the  mass  of  an  O2  molecule,  the  gas  volume,  the  binding  energy  of  O2,  and  e° 
the  nonthermal  potential  energy  of  the  O2  molecule  in  the  gas. 

3.6.  The  Evaluation  of  n+ 


Using  Eq.  (3.3.3),  which  describes  equilibrium  at  the  anode,  along  with  Eq.  (3.3.5)  for  the 
chemical  potentials,  Eq.  (3.3.3)  becomes 


,  9„  9, 

In —  =  In — 


N. 


N 


.  in(-) 


i.^3 
9 


=  In 


9c9g 


(3.6.1) 


or  since  (from  Eq.  (3.3.6)  and  Eq.  (3.4.14))  Na=Nq=l  and  Ng=Wj, 
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= 


w. 


(3.6.2) 


In  an  entirely  similar  manner  for  Eq.  (3.3.4),  one  has 


6  2  3/2  2  6 

In  — -  =  In  - 


2-.  3/2 

w,N^ 


(3.6.3) 


Setting  the  logrithmic  arguments  equal  in  Eq.  (3.6.3)  and  solving  for  N+N®  gives 

2  6>^3/2  2 

6  2  3/2 

^q^c<lx 


(3.6.4) 


The  quantity  q^/Wj  may  be  solved  for  in  Eq.  (3.6.2)  and  substituted  into  Eq.  (3.6.4),  and  since 
charge  neutrality  is  specified,  one  can  set  N  =3N+,  with  the  result 


2  6„ZI2 

XT  n  Q  Q 

3\Ny  = 

2  3/2 

<la^x 


(3.6.5) 


The  partition  functions  given  in  Section  3.5  may  now  be  substituted  into  Eq.  (3.6.5),  after 
which  one  finds 


8  ^  T  r  0  ^ 

X.2  f  e  a  .6^ 


(#t>'^^vexp[-P(6Ve;.2e!.6€“-2e:-1.5e;-1.5e:)] 

X  fl  e 

27t/K  kT 

^  J  Jt^9/4 


(3.6.5) 


In  this,  the  temperatures  T  apply  to  the  components  indicated  by  the  subscripts,  and  q^^  is  defined 
as 


9  ^  2r 

=  %  - 

(=1 

1  -e 


e. 


(3.6.6) 


Equation  (3.6.5)  may  be  divided  by  V*  to  convert  the  total  numbers  N  into  number  densities  n. 
Moreover,  because  of  the  gem-like  quality  of  AI2O3,  the  quantities  that  depend  only  on  the 
volume  of  the  electrolyte  can  be  viewed  as  constant  even  at  detonation  pressures.  These 
quantities  include  the  e’s  associated  with  AI2O3  and  its  enclosed  ions,  electrons,  and  vacancies. 
They  may  be  factored  in  Eq.  (3.6.5)  into  a  constant  along  with  other  constants  such  as  the  binding 
energy  of  O2,  n\,  h,  etc.  Thus,  the  solution  of  Eq.  (3.6.5)  for  n+  is 


•p  -Bp  3/16 ,  3/4  f  P,  1 

n,  =  Vn^e  exp[-^(-2e^-1.5e,)] 


1/8 

9ev 

^9/32 


(3.6.7) 


Here,  p  =  e*+e“+2e“+6e° -1.5e*  and 
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iTlm  k 

p  _  -t  y9B2 


(3.6.8) 


are  constants.  The  quantity  n+  in  Eq.  (3.6.7)  is  the  same  as  n+(/)  in  Eq.  (3.1.11)  and  can  be 
reasonably  evaluated  except  for  the  constant  p,  which  can  be  experimentally  determined. 
Attention  is  turned  next  to  the  cationic  mobility  i|;^. 

3.7.  The  Evaluation  of 

The  cationic  mobility  is  related  to  the  velocity  of  cations  parallel  to  the  walls  of  the 
turmels  described  in  Section  3.4.  The  concepts  and  notation  in  that  development  are  used  here  to 
find  i|;^ 

Consider  a  cation  situated  within  a  tunnel  site.  The  probability  that  it  has  a  kinetic  energy 
parallel  to  the  tunnel  walls  with  some  velocity  component  toward  the  cathode,  which  is  greater 
than  the  barrier  height  et,  is 


/ 


2m 

e  dp 


Pr  " 


_ 


/ 


e  ‘dp 


(3.7.1) 


The  average  value  of  the  velocity,  v ,  in  the  same  direction  for  energies  greater  than  is 


kT 


V  = 


/ 


-p-£l 

^  e  ^"'‘dp 


m  kT 


/ 


e  ‘dp 


(3.7.2) 


Tunnels  are  assumed  to  extend  fi-om  anode  to  cathode,  but  are  not  necessarily  parallel  to 
the  x-axis  of  the  electrolyte.  If  the  tuimel  consists  of  segments  forming  various  angles  at  their 
joints,  the  velocity  in  Eq.  (3.7.2)  must  be  ultimately  corrected  to  give  the  value  in  the  direction  of 
the  current  (parallel  to  the  x  axis).  For  example,  if  lines  parallel  to  all  tunnel  segments  at  some 
value  of  X  form  random  angles  with  the  x-axis,  varying  between  nil  and  -Tr/2,  then  4  v/tu^  will 
give  the  average  velocity  with  which  cations  pass  firom  one  site  to  another  in  the  positive  x- 
direction.  For  the  moment,  assume  this  is  the  case.  Then  the  number  of  cations/sec  per  cm^,  i+, 
passing  through  a  plane  perpendicular  to  the  x-axis  toward  the  cathode  at  position  x  is 


(3.7.3a) 


In  the  product  p^v ,  the  nonintegrable  numerator  in  Eq.  (3.7.1)  and  denominator  in  Eq. 
(3.7.2)  cancel.  The  remaining  integrals  are  easily  evaluated,  giving 
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kT  xl/2  kT 

iTZm 

C 


(3.7.3b) 


If  an  electric  field  e  points  in  the  positive  x-direction,  the  barrier  heights  are  lowered  toward 
positive  X  due  to  the  electric  potential.  The  barrier  to  the  right  (greater  x)  of  the  point  x  in  any 
tunnel  is  lowered  by  3ee(s+Sb)  relative  to  the  barrier  to  the  left.  In  this  case,  using  Eqs.  (3.7.3a) 
and  (3.7.3b)  the  ionic  current  moving  toward  positive  x  is 


Jr 


2'Km 


kT 


(3.7.4) 


The  current  toward  negative  x,  j,,  is  identical,  except  there  is  no  electric  term  in  the  exponential. 
The  net  ionic  current  j+=Jr+Ji  is,  therefore. 


( 


4  ,  kT  -1/2 


3ee(s+s^) 


(3.7.5) 


One  expects  the  electric  potential  to  be  small  relative  to  kT,  allowing  a  two-term  expansion  of 
the  exponential.  This  gives 


J\ 


3en 

C 


4  ^  kT 

7^2  2  Tim 


•) 


1/2 


3ee(s+irj) 

kT 


(3.7.6) 


Because  of  the  definition  of  the  mobility  i|;^  (the  cation  velocity/unit  electric  field),  one  also  has 

]\  =  3en^i|r"e  .  (3.7.7) 

Setting  Eq.  (3.7.6)  and  Eq.  (3.7.7)  equal  and  solving  for  i|;^,  the  result  is 

tj;^  =  ^is+s,)(2%mJcT)-^'h  .  (3  7  8) 

71^  ^ 

This  result,  proportional  to  is  somewhat  different  for  the  usual  form^^  that  is 
proportional  to  T'.  The  change  stems  directly  fi-om  the  fact  that  the  cations  are  allowed  to  move 
freely  between  turmel  barriers,  whereas  in  the  usual  derivation,  they  oscillate  harmonically  in  a 
site  before  jumping  to  the  next. 


One  should  note  that  the  assumed  tunnel  structure  (random  joint  angles,  etc.)  may  not  be 
accurate,  but  differences  in  configurations  are  expected  to  change  i];'^  only  by  a  multiplicative 
constant. 

With  Eq.  (3.7.8)  for  and  Eq.  (3.6.7)  for  n+,  Eq.  (3.1.1 1)  for  K  is  evaluated  except  for 
the  quantity  r4(0),  which  is  assumed  here  to  be  negligably  small.  This  is  because  for  n+(0)  to  be 
significant,  it  would  require  A1  to  be  separated  from  AI2O3  at  the  anode  with  no  oxygen  present. 
This  is  a  highly  improbable  reaction,  so  1^(0)  will  be  neglected.  In  writing  a  final  expression  for 
K,  the  quantity  ng  has  to  be  set  equal  to  1/Qg,  constant  terms  are  again  grouped  for  convenience, 
and  the  electrolyte  temperature  Tg  is  assigned  to  the  expression  for  tj;'^.  The  result  for  K  is 
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where 


Pi 


,5/4 


2€\l.5el 


„  -p  kTf  3/16  q3/4.  ^  e  .1/8  SkT. 


3/A  9132 

^  a  ^  X 


r,  =  e;“it(*+.,)(2ic»,j'”r 

IC 


(3.7.9) 


(3.7.10) 


and 


Pi  =  P+S  •  (3.7.11) 

The  Oj  number  density  n^  is  proportional  to  total  pressure  p.  Note  that  the  3/16  power  is 
close  to  the  1/6  power  heuristically  attached  to  explosive  experiments  by  Guirguis  and  Miller.^* 
The  temperature  structure  of  Eq.  (3.7.9)  is  unusual.  If  all  elements  were  at  the  same  temperature 
(Tf=Te=Ta=Tx),  the  exponential  coefficient  would  vary  as  T™^,  where  Tf  is  the  product  gas 
(fluid)  temperature  in  Section  2.  This  will  be  the  case  in  the  detonation  model  of  Section  2. 
Because  AI2O3  forms  at  the  AljOj/fluid  interface  and  there  dissipates  its  excess  energy,  it  is 
reasonable  to  assume  the  bulk  AI2O3  and  fluid  temperatures  are  roughly  equal.  Moreover,  at  the 
AI/AI2O3  interface,  one  expects  at  least  some  small  depth  of  A1  to  be  at  the  same  temperature  as 
AI2O3 .  Thus,  in  this  model,  Tj=Tj=Tx=Tf. 

Clearly,  the  derivation  of  Eq.  (3.7.9)  neglects  the  intermediate  oxidation  states  of  Al, 
which  must  occur  before  the  formation  of  AI2O3.  These  have  been  assumed  to  be  sufficiently 
short-lived  so  as  to  have  little  effect  on  the  rate  constant.  This  subject  will  be  examined  further 
in  a  future  Phase  II  effort. 


Before  leaving  this  subject,  it  should  be  mentioned  that,  for  convenience,  the  e’s  in  the 
exponents  of  Eq.  (3.7.9)  have  been  expressed  in  terms  different  from  the  heats  of  formation  in 
Section  2.  Nevertheless,  the  two  are  directly  related,  and  the  transformations  are  easily  made  if 
one  wishes  to  do  so.  For  example,  the  heat  of  formation  of  AI2O3,  ^  ,  at  0°K  and  1 

atmosphere  is  per  mole,  ^ 


where  A  is  Avogadro’s  number,  and  the  small  pV  terms  have  been  neglected. 


3.8  The  Al  Reaction  Fraction  A,  As  a  Function  of  Time. 


The  equations  developed  thus  far  for  the  reaction  of  Al  have  been  restricted  to  Al  slabs. 
In  this  calculation,  identical  spheres  of  Al  are  assumed  in  a  mixture  of  Al  and  AP.  Immediately 
behind  the  detonation  front,  these  spheres  will  have  radius  r^,  smaller  than  their  initial  radius  r^. 
In  the  bath  of  hot  product  gas,  AI2O3  will  grow  on  the  sphere.  Let  the  composite  sphere  have 
radius  R;  the  enclosed  Al  have  radius  r^;  r  be  any  specific  radius  in  the  composite;  the  current 
density  in  the  electrolyte  be  j(r). 

j(r)  will  be  dominated  by  the  Al  surface  area  and  decrease  with  increasing  r  such  that  the 
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total  current  is  constant.  That  is. 


Kr)  =  —KR)  . 

Using  Eq.  (3.1.8)  for  j(r)  (with  v|r.»  il;+),  one  has 

dn, 

—  KR)  = 


(3.8.1) 


R^  dr 

Integration  of  Eq.  (3.8.2)  between  r^  and  R  yields,  after  some  rearrangement  (and  neglecting 
ih(ra)), 

3 


(3.8.2) 


KR)  =  \2kT^)f 


nK^) 

R 


a  . -1 


R 


(3.8.3) 


The  time  derivative  of  the  A1  radius  is  (Qa=volume/atom  of  Al), 


dt 


'■a 

=  m 


3e 


R 


Q 

a 

3e 


(3.8.4) 


Or,  using  Eq.  (3.8.3)  for  j(R), 

dt  3e 


R 


R- 


(3.8.5) 


Using  the  formula  for  K  in  Eq.  (3.1.11)  (with  n+(0)=0),  one  can  write  Eq.  (3.8.5)  in  the  form 

2Q 


dr 

a 

dt 


3( - -)K 

Q 


R(i~) 

R^ 


From  Section  2,  the  reacted  fraction  of  Al  is  X=l-r^/ih.  Therefore, 


^  =  -3 

dt 


r^  dr 


.3  dt 


(3.8.6) 


(3.8.7) 


Substitution  of  Eq.  (3.8.6)  into  Eq.  (3.8.7)  gives 

dX 


dt 


=  IH-^)K 


n  2,  2  3 

Q.  r,R 


(3.8.8) 


Now  note  that  the  volume  increase  of  AljOj  is  times  the  volume  change  of  Al. 

This  leads  to 


d3  3  .  3  3. 

R  fn  =  -  (r^-r^)  . 

<7  2Q 


(3.8.9) 
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From  Eq.  (3.8.9)  and  the  definition  of  X,  one  can  show  that 


Q  rl 

e  ^  « 

iQ 

a 


and 


Q 

R  =  r.[l-X+-^X] 
"  2Q 


1/3 


(3.8.10) 


(3.8.11) 


Using  Eq.  (3.8.9)  for  the  described  quantity  in  (3.8.8)  and  replacing  the  resulting  R  with  Eq. 

(3.8.1 1),  one  finally  obtains  Eq.  (3.8.7)  in  the  form  (note  r^/r^  =  (1  -X)^'^) 

2Q  Q 

9( — -fKil  -X+—^Xf^ 

^  ^  (3.8.12) 

* 

With  K  given  by  Eq.  (3.7.9)  and  2QJQ^  by  2V,W2  (from  Section  2),  this  completes  the  rate 
equation  for  the  Al  reaction.  The  reaction  is  assumed  to  occur  only  with  the  Oj  in  the  explosive 
gases.  Reactions  between  Al  and  other  species  are  expected  but  are  considered  much  less 
probable  because  they  are  energetically  less  favorable.  Except  for  the  variation  with  gas 
pressure,  Eq.  (3.8.12)  does  not  resemble  the  heuristic  fit  for  dX/dt  mentioned  earlier.^*  In  the 
latter,  dX/dt  is  found  proportional  to  il-X)\  which  it  clearly  is  not  in  Eq.  (3.8.12).  K,  however, 
has  a  significant  temperature  and  volume  dependence,  so  in  the  end,  Eq.  (3.8.12)  may  show  the 
observed  behavior. 

The  sensitivity  of  the  reaction  rate  with  temperature  may  explain  the  reaction  delays 
directly  and  indirectly  observed  in  many  explosive  experiments.  Immediately  behind  the  shock 
front  of  a  detonation,  the  temperature  may  not  be  high  enough  to  cause  rapid  reaction. 
Nevertheless,  as  reaction  proceeds,  temperatures  may  increase  to  the  point  of  rapid  acceleration 
of  reaction.  This  type  of  process  would  generally  seem  like  a  reaction  delay. 

3.9.  The  Reaction  Rate  of  AP 

Consider  a  system  composed  of  uniform  AP  spheres  of  radius  a.  At  the  front  of  a  plane 
detonation  wave  passing  through  the  system,  the  spheres  will  be  deformed  into  polyhedra  with  no 
intervening  pore  spaces.  Let  these  be  approximated  by  spheres  of  radius  ay.  The  polyhedra 
surfaces,  due  to  friction  generated  at  the  interfaces  and  due  to  pore  collapse,  should  be  at  a  very 
high  temperature  T^.  This  is  assumed  to  cause  immediate  reaction  at  the  surface  with  the 
production  of  products  at  temperature  Tf.  The  unreacted  molecules,  at  temperature  T^p  adjacent 
to  the  products,  will  receive  thermal  energy  from  them,  but  it  will  take  a  certain  critical  time 
for  sufficient  energy  to  reach  the  intramolecular  vibrational  modes  of  the  AP  molecules  to  cause 
reaction  and  produce  equilibrium  products.  This  time  delay  x,.  will  precede  the  transmission  of 
reaction  from  one  molecular  surface  layer  to  the  next.  If  the  molecular  separation  is  A,  then  the 
reaction  proceeds  at  the  rate  dp/dt=-A/x£,  where  p  is  the  radius  of  the  sphericalized  polyhedra  of 
imbumed  AP  and  t  is  time.  The  fraction  of  AP  reacted  at  any  time  is  X,=l-pVaH.  One  can 
combine  these  last  two  equations  to  find,  dX/dt,  and  obtain 
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^  _  3(1 
dt 


(3.9.1) 


The  critical  time  delay  of  devolves  from  three  basic  processes.  The  first  is  a  process  of 
heat  conduction  (thermal  diffusion)  by  which  energy  is  transferred  from  the  product  gas  at  Tf  to 
the  unreacted  AP  at  T^p,  raising  the  temperature  of  the  latter.  The  second  is  a  thermal  relaxation 
process  in  which  energy  is  transferred  from  the  acoustic  vibrational  modes  into  the  vibronic 
(intramolecular)  modes  of  unreacted  AP.  The  third  is  a  process  by  which  the  AP  is  reduced  to 
products  once  reaction  has  begun.  Each  of  these  processes  takes  time  and  contributes  to  x^.  The 
time  Xd  for  the  first  process  is  approximately  PA^  where  P  is  the  reciprocal  of  the  diffusivity  of 
AP.  At  standard  temperature  and  pressure  (STP),  P~601  sec/cm^  and  A«5.8A“,  which  makes  Xp 
of  the  order  of  2  ps.  This  time  is  too  small  to  govern  the  reaction  rate  of  AP,  which  in  9p-sized 
particles  is  known*  to  react  in  about  0.5  psec  under  detonation  conditions.  If  x^  were  equal  to  x^, 
this  reaction  would  occur  in  15  ns.  One  should  therefore  expect  that  the  unreacted  AP  adjacent 
to  product  gas  reaches  highly  elevated  temperature  well  before  reaction  proceeds  to  products. 


The  second  process  is  governed  by  a  relation  of  the  type'** 


T  -T 

•'  c  ^  AP 


^-202y, 

^  02 


(3.9.2) 


where  T^,  is  some  critical  temperature  associated  with  the  energy  an  AP  molecule  must  reach 
before  reaction  begins  promptly  (in  say  less  that  1  ps).  0  is  the  Debye  temperature  of  the  acoustic 
modes  in  AP,  6  is  the  ratio  of  the  Einstein  temperature  of  the  lowest  frequency  vibronic  mode, 

0E,  to  0(i.e.,  6=0e/0),  and  ^  is  a  constant  characteristic  of  AP. 


When  both  sides  of  Eq.  (3.9.2)  are  multiplied  by  k,  and  Xr  is  replaced  by  t,  it  represents 
the  energy  poured  into  the  lowest  frequency  vibronic  mode  in  time  t  by  a  neighboring  molecule  at 
temperature  Tf.  The  process  stops  when  the  energies  are  in  balance.  One  expects  this  relatively 
slow  process  to  determine  the  rate  at  which  energy  transfers  into  vibronic  modes. 


For  6!  in  Eq.  (3.9.2),  one  may  use  the  Stirling  formula,  which  is  reasonably  accurate  even 
for  6=1.  So  set 

6!  =  6®e  ■®(27i:6)*'2  (3  9  3^ 


The  temperature  0e  in  delta  may  be  determined  from  spectral  data  on  AP,  and  0,  which  varies 
with  volume,  can  be  expressed  in  terms  of  the  Gruneisen  parameter  Ya  for  the  highest  frequency 
AP  acoustic  modes.  By  definition. 


-41n0) 

d{\nv^p) 


(3.9.4) 


When  integrated,  this  gives 

Y 

0  =  00  exp[-  f  dv^p]  . 

V  V.p 


(3.9.5) 
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Reasonable  expressions'”’'*®  for  may  be  derived  from  the  equation  of  state  of  AP."  Thus,  6  and 
0  can  be  determined  as  a  function  of  v^p. 

Solving  for  Tr  in  Eq.  (3.9.2),  using  Eq.  (3.9.3)  for  6!,  one  has 

r  1/2 

Tr  =  [6^e  ■®(27t6)”2  '^c 

02 

This  represents  the  time  required  for  a  molecule  of  AP  at  T=Tap  to  react  when  adjacent  to 
molecules  at  T=Tf.  The  quantities  ^  and  T^  are  constants  related  to  AP  only  and  must  be 
experimentally  determined.  It  should  be  noted  that  the  relation  Eq.  (3.9.2)  is  classically 
developed.'**  Except  for  a  multiplicative  constant  (near  1),  the  quantum  mechanical  analogue  is 
identical. 

The  third  process  contributing  to  involves  the  time  Xg  necessary  for  reaction  to 
proceed  to  equilibrium  products  at  Tf,  when  the  average  molecular  temperature  is  sufficient  to 
initiate  reaction  at  the  begiiming  of  the  process.  'This  time  should  be  some  multiple  of  the  time  it 
takes  a  molecular  fragment  to  collide  with  another  molecule.  In  AP  at  STP,  the  separation 
between  molecules  is  about  5. 8 A,  but  the  molecules  themselves  have  dimensions  (sphericalized 
diameters)  of  about  3.36  A.  So  the  average  distance  to  collision  for  a  fragment  is  about  2.5A  at  a 
velocity  of  roughly  (kT/mf)*'",  where  mj  is  the  fragment  mass. 

Taking  Tf  to  be  (a  conservative)  1000°K,  the  time  between  collisions  for,  say,  an  NH4 
fragment  is  about  0.4  ps.  If  ten  collisions  are  required  to  convert  the  fragment  components  into 
parts  of  equilibrium  products,  then  only  4  ps  is  the  requisite  time.  In  view  of  this,  it  seems 
reasonable  to  assume  that  the  first  and  third  processes  contributing  to  Xg  are  small  compared  to 
the  second,  which  is  in  the  approximate  range  10-120  ps^*  '*®  and  that  Xg  can  be  represented  by 

\  ’  (3.9.7) 

where  x'-  Xo+Xg  and  is  treated  as  a  constant  with  a  value  between  about  2  ps  and  10  ps. 
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4.  CONCLUSIONS  AND  SUMMARY 

4.1.  Technical  Results 

Theoretical  work  was  performed  to  formulate  a  first  generation  model  for  nonideal 
detonation  that  is  based  on  fundamental  physical  principles  and  can  be  used  to  investigate  the 
nonideal  detonation  process  in  explosive  compositions  containing  Al. 

In  the  portion  of  this  work  presented  in  Section  2,  we  formulated  the  thermo¬ 
hydrodynamics  that  is  needed  to  treat  reacting  AP/Al  and  CHNO/Al  compositions  and  account 
for  the  decomposition  of  these  explosives  and  the  subsequent  reactions  of  their  products  with  Al. 
We  also  considered  the  relationship  between  the  steady-state  reaction  zone  length  and  reaction 
time  in  AP,  and  developed  a  prototype  model  for  steady-state,  nonideal  detonation  in  a  CHNO/Al 
composition. 

In  our  Phase  I  approach  to  reactions  in  these  AP/Al  and  CHNO/Al  compositions,  single¬ 
reaction  coordinates  were  used  to  describe  the  decompositions  of  the  AP  and  the  CHNO 
explosive,  and  a  single  reaction  coordinate  was  used  to  describe  the  reactions  of  their  products 
with  Al  to  form  either  a^Al203  (corundum)  or  amporphous  y-AljOj.  The  second  reaction 
coordinate  for  the  AP/Al  composition  was  assigned  to  the  reaction  of  Al  with  Oj;  the  second 
reaction  coordinate  for  the  CHNO/Al  composition  was  assigned  to  the  reaction  of  Al  with  HjO 
and  C,  with  the  Hj  produced  from  the  H2O  reacting  with  the  C  to  form  C2Hg.  Although  these 
reaction  schemes  are  clearly  limited,  they  are  more  realistic  than  some  empirical  treatments  of  Al 
reactions  in  hydrocodes,  because  they  account  explicitly  for  the  conservation  of  mass  as  the  Al 
removes  oxygen  from  the  detonation  products. 

The  mass-balance  equations  for  these  reaction  schemes  were  used  to  construet 
constitutive  equations  for  the  specific  internal  energies  of  our  reacting  AP/Al  and  CHNO/Al 
mixtures,  with  each  of  their  components  governed  by  its  own  equation  of  state  (EOS). 
Generalized  Mie-Gruneisen  equations  of  state  were  used  to  describe  the  condensed  components, 
and  the  fluid  detonation  products  were  assumed  to  be  polytropic  for  convenience.  Because  heat 
conduction  was  not  treated  explicitly  in  constructing  these  constitutive  relationships,  assumptions 
were  required  for  their  closure.  Assumptions  were  made,  (1)  that  no  appreciable  amount  of  heat 
is  transferred  into  the  interior  of  explosive  and  Al  grains,  and  (2)  that  the  heat  transfer  process 
between  the  AljOj  covering  the  Al  grains  and  the  detonation  products  is  efficient.  As  a  result  of 
the  first  assumption,  burning  explosive  and  reacting  Al  grains  in  an  expansive  flow  behind  the 
shoek  are  constrained  to  their  own  isentropes  and  can  adequately  be  described  by  their  respective 
energy,  pressure,  volume  equations  of  state.  As  a  result  of  the  second  assumption,  the  AI2O3  and 
detonation  products  attain  thermal  equilibrium  and  the  temperature  equations  of  state  of  these 
components  are  required  to  formulate  this  condition. 

Equations  of  state  for  the  components  in  our  reacting  AP/Al  composition  were 
constructed  to  make  our  constitutive  equation  for  this  system  praetical.  Our  equation  of  state  for 
Al  was  based  on  the  Morse  potential,  and  with  values  of  the  latent  heat  of  vaporization,  the 
specific  heat  at  constant  volume,  and  the  Gruneisen  parameter  for  Al  from  the  literature,  predicts 
shock  wave  states  that  are  in  excellent  agreement  with  those  determined  experimentally  up  to 
400  kbar.  No  equation  of  state  was  constructed  for  liquid  Al,  however,  because  melting  of 
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condensed  components  was  excluded  from  our  first-generation  model. 

The  equation  of  state  for  AP  was  based  on  shock  wave  and  velocity  of  sound  data  and  the 
assumption  that  its  specific  heat  at  constant  volume  and  its  Gruneisen  parameter  are  constants. 
The  unexpected  high  values  for  shock  temperature  calculated  in  the  300-400  kbar  region, 
however,  suggests  that  the  validity  of  these  assumptions  must  be  considered  in  any  future  work 
on  AP/Al  compositions. 

The  first  generation  EOS,  constructed  for  a-AljOj  and  amorphous  y-AljOj,  were  based 
on  shock-wave  data  for  corundum  and  the  assumptions  that  these  oxides  behave  as  linear  elastic 
materials  with  constant  specific  heats  at  constant  volume  and  constant  Gruneisen  parameters  in 
the  pressure  range  of  interest.  Literature  values  of  the  density,  the  standard  heat  of  formation  and 
standard  entropy  for  a-Al203,  and  of  the  densities  and  standard  heats  of  formation  of  amorphous 
y-AljOj  were  used  to  evaluate  Gibbs  free  energies  at  different  pressures,  and  thereby  obtain  an 
estimate  of  the  standard  entropy  of  amorphous  Y-AI2O3  and  formulate  a  criterion  for  the 
formation  of  either  a-Al203  or  amorphous  Y-AI2O3  in  our  reacting  AP/Al  mixture.  It  is  important 
to  note  here  that  this  description  of  AI2O3  is  preliminary  and  needs  to  be  improved  in  any  future 
work  that  includes  melting. 

Our  construction  of  a  temperature  equation  of  state  for  the  polytropic  fluid  component 
was  based  on  a  complete  equation  of  state,  which  provides  the  specific  energy  relationship 
already  used  for  the  detonation  products,  the  thermodynamic  identities  for  temperature  and 
pressure,  and  the  assumption  that  the  specific  heat  at  constant  volume  is  constant.  Here  again  the 
description  of  the  fluid  component  must  be  improved  in  any  future  work  on  AP/Al  and 
CHNO/Al  compositions. 

Certain  aspects  of  steady-state  detonation  were  considered  to  provide  a  better 
understanding  of  the  dependence  of  the  nonideal  detonation  process  on  reaction  rates.  Equations 
for  estimating  upper  and  lower  limits  for  a  steady-state  reaction  zone  length  were  derived;  shock 
wave  data  and  the  results  of  a  TIGER  code  calculation  for  AP  were  used  to  obtain  equations  for 
the  upper  and  lower  limits  of  its  reaction  zone  length  in  terms  of  its  reaction  zone  time.  Our 
prototype  model  for  steady-state  nonideal  detonation  was  obtained  by  subjecting  our  previous 
treatment  of  the  CHNO/Al  composition  to  simplifying  assumptions  about  the  condensed 
components.  The  equations,  derived  for  Chapman-Jouguet  parameters  in  terms  of  the  reaction 
coordinate  describing  the  reaction  of  A1  with  detonation  products,  demonstrate  clearly  the 
dependence  of  steady-state  nonideal  detonation  on  this  reaction  and  the  reason  for  developing  a 
realistic  expression  for  its  rate. 

In  the  portion  of  our  theoretical  work  presented  in  Section  3,  we  formulated  rate 
expressions  for  the  oxidation  of  A1  and  the  decomposition  of  AP,  which  are  needed  to  model 
nonideal  detonation  in  our  AP/Al  compositions. 

For  the  derivation  of  the  rate  equations  for  the  A1  reaction ,  the  AI/AI2O3/O2  system  is 
treacted  as  a  galvanic  cell  in  which  A1  ions  at  the  AI/AI2O3  interface  are  injected  into  the  AI2O3 
electrolyte  and  migrate  toward  the  AI2O3/O2  interface  where  AI2O3  is  created. 

In  the  construction  of  the  rate  model,  it  is  assume  that  at  the  high  pressures  and 
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temperatures  attained  during  the  detonating  process,  P-AljOj-like  materials  are  present,  which 
are  known  to  exhibit  anomalously  high  ionic  conductivity  because  of  the  formation  of  tunnels 
through  which  ions  can  more  easily  pass.  The  equations  describing  the  reaction  process  are 
presented  and  the  rate  constant,  ionic  conductivity,  and  mobility  are  calculated.  This  results  in  a 
rate  constant  that  is  an  Arrhenius  term  multiplied  by  a  factor  proportional  to  T™p^^’®. 

It  is  concluded  that 

a.  AI2O3  forms  primarily  at  the  AI2O3/O2  interface.  This  is  based  on  the  excessive  time  it 
would  take  for  the  heat  of  reaction  to  diffuse  from  the  AI/AI2O3  interface  into  the  fluid  products 
of  the  AP  reaction. 

b.  Tunnels  as  ion  conduits  (or  some  similar  scheme)  are  necessary  in  order  to  account  for 
the  high  conductivity  required  to  produce  the  observed  reaction  rates. 

The  reaction  rate  of  AP  is  thought  to  be  dominated  by  thermal  relaxation  processes,  the 
volume  and  temperature  dependences  of  which  are  estimated.  The  time  it  takes  for  energy  to 
transfer  from  the  acoustic  vibrations  of  the  AP  molecule  into  its  vibronic  modes  is  seen  to  be 
long  compared  to  other  processes.  It  is  assumed  that  reaction  begins  with  the  excitation  of 
vibronic  modes  and  that  the  thermal  relaxation  time  is  governed  by  energy  transfer  from  acoustic 
modes  into  the  lowest  frequency  vibronic  (“doorway”)  mode  from  which  energy  is  fed  much 
more  quickly  into  other  vibronic  modes. 

The  time  necessary  to  excite  a  critical  vibronic  mode  therefore  is  inversely  proportional  to 
the  rate  at  which  energy  enters  the  doorway  mode.  This  rate  is  dependent  on  the  ratio  of  the 
doorway  frequency  to  the  highest  frequency  of  the  acoustic  modes  and  is  a  minimum  if  the  ratio 
is  one.  The  rate  is  also  volume-dependent  because  of  the  change  of  the  acoustic  mode 
frequencies  with  volume. 

The  rate  equations  developed  here  for  the  production  of  both  AI2O3  and  the  explosive 
products  of  AP  are  not  especially  specific  to  these  substances.  They  should  apply  as  well  to  a 
number  of  explosives,  oxidizers,  and  metals  when  appropriate  substitutions  for  the  values  of  the 
key  parameters  are  made. 

4.2.  Summary  of  Progress  vs.  Proposed  Tasks 

The  work  described  in  Sections  2  and  3  of  this  report  followed  in  general  along  the  lines 
described  in  the  Phase  I  Scope  of  Work  in  our  proposal  of  14  July  1994,  but  deviated  in  some 
important  ways. 

Task  1  of  the  Scope  of  Work  sought  to  “develop  a  preliminary  model  for  transport  and 
kinetics  mechanism[s]  that  determine  the  energy  conversion  and  energy  release  rates  in  a  totally 
nonideal  explosive,  such  as  an  AP/Al  mixture, ...”  Section  3  of  this  report  describes  a  theoretical 
model  of  the  kinetics  of  the  AP,  Al,  and  AI2O3  interactions,  which  is  far  more  complex  and 
sophisticated,  albeit  preliminary,  than  attempted  previously.  The  heat  conduction  processes  are 
handled  implicity  in  both  Sections  2  and  3. 
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Task  3  sought  to  “examine  different  equations  of  state  models  as  candidates  for 
representing  the  products  of  metal/fuel/oxidizer/CHNOF  explosive  mixtures  and  metal 
fuel/oxidizer  mixtures  ...”  Inasmuch  as  only  an  AP/Al  mixture  was  to  be  considered  (by  Task  1), 
the  most  important  detonation  products  are  N2,  Oj,  H2O,  CI2,  and  a-  and  amorphous  AI2O3 
(solid).  Section  3  contains  the  equations  of  state,  both  E,P,V  and  T,P,V,  derived  for  AI2O3. 
Equations  of  state  for  A1  and  AP  as  well  as  polytropic  equations  of  state  for  CHNO  detonation 
products  that  react  with  A1  are  also  derived. 

Task  2  was  to  “insert  the  preliminary  model  of  Task  1  into  a  hydrodynamic  code  and 
calculate  the  failure  diameter  of  a  charge...”;  Task  4  was  “use  a  1-D  hydrocode  to  calculate  the 
underwater  bubble  expansion  to  the  first  maximum... using  the  results  from  Tasks  1  and  3  ...” 
Early  on,  we  started  work  on  Task  4  but  quickly  dropped  it,  as  well  as  Task  2.  We  decided  to 
dedicate  almost  the  entire  effort  to  the  much  more  important  problems  of  developing  the  physical 
models  of  the  kinetics  of  the  AP/Al  reaction,  the  development  of  the  needed  temperature 
equations  of  state  for  reactants  and  products,  particularly  to  accomodate  our  assumption  that  the 
AI2O3  and  detonation  products  are  in  thermal  equilibrium,  and  to  describe  a  phenomenological 
model  based  on  fundamental  physical  principles. 

The  insertion  of  the  physical  models  developed  in  Phase  I  into  hydrocodes  has  been 
deferred  to  Phase  II  because  of  time  constraints. 

4.3.  Recommendations. 

The  thermo-hydrodynamic  model  presented  in  this  report  is  based  on  fundamental 
physical  principles  and  provides  a  firm  foundation  for  a  first-generation  hydrocode  for 
investigating  the  nonideal  detonation  process  in  explosives  containing  Al.  In  addition,  it  is  more 
realistic  than  many  currently  used  in  hydrocodes  because,  (1)  untenable  assumptions,  such  as  the 
attainment  of  thermal  equilibrium  in  the  reacting  explosive  mixture,  have  been  removed,  and  (2) 
the  need  for  empirical  fits  of  the  rate  equations  has  been  eliminated.  Nevertheless,  our  thermo¬ 
hydrodynamic  model  is  clearly  limited  by  restrictive  assumptions,  imposed  by  Phase  I  time 
constraints. 

Proposed  tasks  for  a  Phase  II  program,  which  follow  naturally  from  the  work  performed 
in  Phase  I  and  are  needed  to  remove  the  latter’s  restrictive  assumptions,  are  listed  below: 

1.  Incorporate  the  first-generation  model  into  a  hydrocode,  and  validate  it  by  performing  a 
series  of  1 -dimensional,  hydrodynamic  calculations  for  AP/Al  mixtures  and  comparing  their 
detonation  velocities  with  experimental  data. 

2.  Extend  our  treatment  of  AP/Al  and  CHNO/Al  compositions  by  including  porosity. 

3.  In  parallel  with  other  recommendations  above  and  below,  incorporate  as  modules  into 
an  existing  2-dimensional  flow  hydrocode,  e.g.,  CALE,  the  first-  and  subsequent-generation 
models,  perform  calculations  for  AP/Al  and  AP/CHNO/Al  mixtures,  and  compare  their  results 
with  experimental  data  available  from  Moby-Dick  experiments  at  NSWC  and  underwater 
detonation  experiments,  including  the  bubble  oscillation  phase. 
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4.  Extend  the  reaction  schemes  developed  for  AP/Al  and  CHNO/Al  compositions  in 
Phase  I,  by  including  the  reaction  of  A1  with  HjO  in  the  AP/Al  system  and  the  reaction  of  A1  with 
COj  in  the  CHNO/Al  system. 

5.  Examine  the  mechanism  of  heat  transfer  in  reacting  explosive  mixtures  containing  Al, 
with  particular  emphasis  on  late  time. 

6.  Remove  the  restrictive  assumptions  that  the  values  of  and  F  are  constants  for 
unreacted  explosives  and  then  construct  more  realistic  EOS  for  unreacted  AP  and  a  prototype 
CHNO  explosive. 

7.  Improve  the  equations  of  state  of  a-Al203  and  amorphous  Y-AI2O3  and  formulate, 
based  on  Gibbs  free  energy  considerations,  a  more  basic  criterion  for  the  conversion  of  one  into 
the  other. 

8.  Consider  melting  of  Al  and  AI2O3  in  the  reacting  AP/Al  and  CHNO/Al  mixtures. 

9.  Replace  the  polytropic  EOS  for  the  fluid  detonation  products  by  a  more  fundamental 
EOS  based  on  statistical  mechanics  of  the  type  developed  and  implemented  by  F.H.  Ree,  W.B. 
Brown  and  H.D.  Jones  and  F.  Zerilli. 

10.  Extend  our  thermo-hydrodynamic  modeling  to  reacting  AP/CHNO  explosive/Al 
mixtures. 

1 1 .  Investigate  the  intermediate  oxidation  states  of  Al,  e.g.,  AI2O2  and  AIO2,  to  ascertain 
if  other  Al  oxides  are  formed  in  rate  controlling  steps  in  the  formation  of  AI2O3. 

12.  Include  secondary  oxygen  sources,  such  as  H2O,  and  CO2,  in  the  reaction  rate 
expression  for  Al,  and  consider  its  oxidation  by  chlorine. 

13.  Examine  the  effect  of  replacing  Al  with  another  metal. 
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APPENDIX  A.  A  MODEL  FOR  STEADY-STATE  REACTION  ZONE  (SSRZ)  TIMES 

Our  treatment  of  steady-state  (SS)  detonation  here  is  based  on  the  hypothesis  that  the 
global  decomposition  of  a  conventional  condensed  explosive  is  pressure  dependent.  In  this  case, 
it  is  significant  to  consider  how^  SSRZ  times  are  governed  by  an  explosive’s  pressure-dependent 
burning  rate  because  the  classical  treatment  of  the  dependence  of  the  SSRZ  on  particle  size  is 
based  on  the  assumption  that  the  explosive  bums  at  a  constant  rate.  Accordingly,  we  develop  a 
prototype  model  for  a  SSRZ  in  order  to  investigate  the  relationship  between  an  explosive’s  SSRZ 
time  and  its  burning  rate.  To  provide  the  necessary  backgroxmd  for  this  development,  we  first 
present  the  conservation  equations  and  Chapman- Jouguet  (CJ)  conditions  for  a  SS  detonation 
wave. 

Equations  for  a  SSRZ 

We  use  the  subscripts  x  and  p  to  denote  the  explosive  and  its  products,  the  subscript  H  to 
denote  the  shocked  state,  and  either  a  superscript  or  a  subscript  0  to  denote  the  standard  state  of 
the  explosive  and  the  detonation  products,  i.e.,  at  pressure  po  and  temperatures  Tq.  Recalling  that 
the  Rankine-Hugoniot  (RH)  jump  conditions  are  first  integrals  of  the  differential  equations 
expressing  the  conservation  laws  when  the  flow  has  attained  a  SS,  we  write  the  equation 
governing  our  SSRZ  as 


p(D-u)=poD, 

(A.1) 

p=poD^(l-v/vo), 

(A.2) 

e-el=pvo(l-v/vQ)/2, 

(A.3) 

with  D  and  u  used,  respectively,  to  denote  the  detonation  velocity  and  the  particle  velocity,  and 
Po~0.  Equation  (A.l)  expresses  the  conservation  of  mass,  Eq.  (A.2)  expresses  the  conservation 
of  mass  and  momentum,  Eq.  (A.3)  expresses  the  conservation  of  mass,  momentum,  and  energy, 
and  the  equations  governing  the  shock  at  the  front  of  the  SSRZ  are  obtained  by  putting  a 
subscript  // on  the  variables  p,u,p,v  and  e  in  Equations  (A.l),  (A.2),  and  (A.3). 

To  continue  this  treatment  of  our  SSRZ,  we  introduce  the  sound  speed  c,  define  the 
adiabatic  index  k  by  the  equation,  c^=kpv,  and  recall  that  the  SSRZ  terminates  at  the  CJ  state 
where  the  flows  become  sonic.  Using  the  subscript  j  to  denote  the  CJ  state,  this  sonic  condition 
is  expressed  by  the  equations 

D-Uj=Cj.  (A.4) 

The  combination  of  the  equations  pfj=pJD,  c]=kjPjVj,  and  Eq.  (A.2)  written  for  the  CJ  state  then 
gives  the  well  known  equations  for  v,u,c,  and  p  in  the  CJ  state  as 


(A.5) 

(A.6) 


Vj/vo=k/(kj+l), 

kjUj=Cj=Dk/(kj+l), 
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Pj=PoD^/(kj+l).  (A.7) 

We  finally  note  that  the  relationship  between  D  and  the  specific  chemical  energy  q 
liberated  by  a  particle  passing  through  the  SSRZ  is  derived  by  combining  the  e=e(p,v)  equation 
of  state  for  the  detonation  products  with  the  following  equation: 

ej-el=D^/2(kj+iy  (A.  8) 

obtained  by  combining  Eqs.  (A.5)  and  (A.7)  with  Eq.  (A.3)  written  for  the  CJ  state. 

Calculations  of  SSRZ  Times 


For  the  sake  of  tractability,  we  assumed  that  the  shocked  explosive  is  incompressible  with 
a  volume  and  a  Gruneisen  parameter  r=v(dp/de).  In  this  case,  with  A,  used  to  denote  the 
mass  fraction  of  decomposed  explosive,  we  can  write  the  mixture  rules  for  the  specific  volume 
and  specific  energy  of  the  explosive  and  its  products  as 

v=(l-A,)v+A,Vp  (A.9) 

e=(l-X)e^+Xep  (A.  10) 


With  the  added  assumption  that Px=Pp-p,  the  e^=e/p,v)  and  ep=ep(p,Vp}  equations  of  state  are 
then  required  to  derive  a  constitutive  equation  for  e  in  terms  of p,v,  and  v.  We  write  the  equation 
for  e^  as 


e 


X 


2 

[i-(i+f^)— ] 
2  r  Vp 


(A.11) 


We  assume,  again  for  tractability,  that  the  explosive’s  products  are  polytropic  with  an  index  y 
and  write  the  equation  for  as 


e 

p 


+ 


Eh. 

Y-l 


(A.  12) 


The  combination  of  Eqs.  (A.9)-(A.12)  allows  us  to  write  the  constitutive  equation  for  e  for  our 
decomposing  explosive  as 


e 


+ 


jpy_ 

Y-l 


+ 


(l-^)Vo[ 


pV 

Y-l 


(A.13) 


with  the  specific  heat  of  decomposition  q  =  -(e^  ~e°),  E=(( y  - 1  )/r  - 1  )v/Vq,  and 
(j)=i  -{i+2fr)Vf/vQ. 


It  is  now  convenient  to  derive  the  relationships  between  k  and  y  and  £F  and  q.  We  first 
use  the  thermodynamic  identity  for  the  sound  speed  and  c^=lq}v,  which  yield 
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de 

dv 


de 


dp 


(A.14) 


and  the  equations  for  the  partial  derivatives  of  e  derived  from  Eq.  (A.  13)  to  obtain  the  equation 
relating  k  and  y  as. 


k  = 


yv 

v+(1-A)VoF 


(A.15) 


Setting  k=kj  and  v=Vj  in  Eq.  (A.15)  when  A=l,  then  shows  that  kj=y  at  the  end  of  the  reaction 
zone.  Combining  the  equations  obtained  by  setting  kj=y  in  Eq.  (A.8)  and  setting  A=l,  e=ep 
v=Vj=yvo  /(y+1),  and Pj=pgD^/(y+\)  in  Eq.  (A.13)  gives  the  well-known  equation  relating  D  and 

q'- 


-  2{f-\)q  .  (A.16) 

The  equations,  relating  X  and  v  and  X  and p  in  the  SSRZ,  can  now  be  derived  from  Eqs. 
(A.2),  (A.3),  (A.13)  and  (A.16)  as  follows.  We  substitute  the  expressions  for  p  and  q  from  Eqs. 
(A.2)  and  (A.16)  into  the  equation  obtained  by  eliminating  e-e®  between  Eqs.  (A.3)  and  (A.13)  to 
derive,  after  some  manipulation,  the  following  equation  relating  X  and  v  in  the  SSRZ: 


[Y-(Y+l)-f 
_ ^0 

1  +2(y  +  l)K(l  -iL)+(y2-l)(|)(l 

Vn  V„ 


(A.  17) 


The  more  interesting  equation  for  our  purpose  here,  that  relates  X  and p  in  the  SSRZ,  is  obtained 
by  combining  Eqs.  (A.  17)  and  (A.2)  as 


- _  =  F(^  ,  ^^1)  . 

l+2F^+(y-l)(j)— 

Pj  Pj 


(A.  18) 


Because  Eq.  (A.  18)  will  be  used  to  calculate  SSRZ  times,  it  is  convenient  to  consider  a 
typical  condensed  explosive  so  that  different  values  of  F  can  be  used  to  evaluate  the  parameters 
V and  4).  We  consider  a  typical  explosive  with  the  parameters:  po=l  .76  g/cc,  D=1A  mm/psec, 
y=3,/7y  =  240.9  kbar,/>H=300  kbar  and  Vh/vo=0.6887.  The  two  cases,  when  F=0  and  (j)=0,  are  of 
particular  interest  here.  When  r=0  and  r=y- 1=2,  Eq.  (A.  18)  becomes 


(^-1)^ 

-  (A.18.1) 

Pj 

with  a=0.7548,  and  when,  (J)=0  and  r=4.425,  Eq.  (A.  18)  becomes 
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1-A  =  ^ 


(A.  18.2) 


Pj 


Before  we  can  derive  an  equation  for  SSRZ  times,  we  must  specify  an  expression  for  the  time 
rate  of  decomposition  of  the  explosive,  dA/dt,  in  terms  of  the  state  variables,  where  t  is  the  time. 
Here,  we  assume  that  our  shocked,  incompressible  explosive  grains  are  spherical,  have  an  initial 
radius  R„  and  bum  with  a  pressure-dependent  rate  governed  by  the  equation 

dR  „ 

(A.19) 

with  n  a  variable  parameter.  In  this  case,  the  equation  relating  A  to  the  grain  radius  R  is, 

A  =  1-(A)3  , 

R. 

I 

and  the  equation  for  the  rate  of  decomposition  follows  by  differentiating  Eq.  (A.20)  as 

dX  ^  _A(A)2^ 

dt  R.  R.  dt 

The  combination  of  Eqs.  (A.21),  (A.20),  and  (A.19)  then  gives  the  equation  for  the  rate  of 
decomposition  of  our  explosive  in  terms  of  A  and  p  as 


(A.20) 


(A.21) 


—  =  —{i-xyAv" . 

dt  R. 


(A.22) 


We  now  use  Eq.  (A.19)  to  define  the  characteristic  time  ttn  as 

_L  =  =  -— (— ) 

a  R.  R^  dt^' 

n  i  i 


(A.23) 


and  rewrite  Eq.  (A.22)  as 


^  _  3(1  -X)^  Pj  p  „ 
dt  Ph  Pj 


(A.24) 


Our  general  procedure  for  calculating  SSRZ  times  is  to  integrate  the  expression  for  dt/dp 
in  the  SSRZ,  obtained  in  terms  of p  from  Eq.  (A.  18)  and  (A.23),  from p/pj~Pf/pj  to  p/pj=\.  Here, 
we  will  exemplify  this  procedure  for  the  case  when  V=0.  When  F=0,  we  equate  the  expression 
for  d/l/dt,  derived  by  differentiating  Eq.  (A.  18.1),  and  the  expression  for  dA^dt,  derived  by 
combining  Eq.  (A.  18.1)  and  (A.23)  to  eliminate  (l-A)‘^^  to  obtain  the  following  equation  for 
dp/dt  in  the  SSRZ: 
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r,  I  p.  dt 

Pi 


a  p: 

«  w 


(A.25) 


Formally  integrating  Eq.  (A.25),  with  P=p/pj  and  Ph=Pi/Pj,  gives  the  equation  for  the  SSRZ  time 


t„  as 


K 

a 


-2P 


n  1 


dP 


-2P 


H 


T'"- 


where  /„  is  the  integral  in  Eq.  (A.26).  Combining  Eq.  (A.23)  and  (A.26)  gives  the  equation 


(A.26) 


-  - 


2R. 


1  t. 


(-)  , 


3(p)\\-aP„y 


(A.27) 


which  allows  us  to  calculate  A„  when  R,  and  n  are  known  and  the  SSRZ  parameters  Ph,  Pj,  and  t„ 
are  also  known. 


We  now  set  n=l  and  evaluate  I,  to  determine  the  relationship  between  tj  and  a,.  To 
reduce  Ij  to  a  more  recognizable  form,  we  change  the  variable  from  P  toy  by  setting y"=P-l,  and 
then  set  w=3  to  obtain  the  equation 


0 


(A.28) 


where  y/,=(P;^l)''^=0.6258.  We  now  use 

/- 

^  1 


y  ^y.  = 

H-y*  ' 


(A.29) 


to  evaluate  this  integral  between  y^O.625 8  andy=0,  and  find  the  following  equation  relating  tj 
and  a,:  t,=\.\?)laj. 


We  now  set  n=2  and  evaluate  the  integral 


(A.30) 


to  determine  the  relationship  between  and  0.2-  An  operation  using  integration  by  parts, 
separation  by  partial  fi-actions,  and  the  change  of  variable  y=P-l  allows  us  to  write  an  equation 
for  I 2  as 


Ph  3  ■ 


(A.31) 
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Equations  (A.26)  and  (A.3 1)  then  give  the  following  equation  relating  and  ttj:  t2=\ 30a!2. 

It  is  clear  that  the  results  of  these  calculations  are  limited  by  our  simple  description  of  the 
explosive  and  its  products.  Nevertheless,  they  are  significant  for  an  explosive  with  a  pressure 
dependent,  decomposition  rate  because  they  show  that  the  parameter  An  in  its  rate  law  can,  in 
principle,  be  determined  when  n,  the  initial  grain  size,  the  ZND  spike  pressure,  and  the  SSRZ 
time  are  known. 
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